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内容简介 

本 -15 介绍优选学方面的一些常觅的方法.全书共分三部分.第一 
部分单囚素 优选法 ，主要是介绍货金分割法、分数法和抛物线法.第二 
部分多 因素优 选法，主要是介绍双因索优选法、最陡上升法，切块法、二 
次回归法 Iri 抛物休法. 第三邱 分附录，介绍其它一 些方法 ， 

本书可供卨等唭校数学系师生以及运筹学工作者参考使用， 
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序 


言 


这本书的原稿完成于六十年代末期，作为向国庆二十周 
年的献礼，但一直未能出版，经过了十年，终于在国庆三十周 
年问世了. 

优选法的试点、推广工作，曾经在毛主席和周总理的关怀 
下，在全国各地广泛地开展起来. 

优选法是尽可能少做试验、尽快地找到生产最优方案的 
方法，也就是现代科学试验的方法.过去总是夸耀某人花了 
几十年功夫，做了几百次、甚至几千次的试验，获得了某一个 
科学成果.自从摆脱了以往排冽组合做试验而出现了优选法 
之后，我们就要反问一句，是否有必要做那么多次试验=花那 
么长的时间才能获得某一成果，当然，花几十年搞试验的精 
神和毅力是令人钦佩的.可是 s 用优选法搞几十年就可以得 
到更多更好的成果. 

优选法首先应用在化工、电子工业，然后逐步推广到各行 
各业.连炸油条省油，做豆腐省豆都用上了.在医学上3卜科 
医生用了优选法能够在病人肠子里找到出血点. 

为了使广大群众能够掌握和运用这个方法，曾经编写了 
几本通俗易懂的小册子.例如， 《优选 法平话》(署名齐念一, 
科学出版社，1971年)、《优选法平话及其补充》（国防工业出 
版社， 1 W 1 年）、《优选法话本》(辽宁人民出版社，年）. 

近年来，我和各省市组织的优选法推广队伍跑遍了祖国 
二十二个省、市、自治区，到一百多个大小城镇，向广大群众介 
绍这个方法，并且推广各地组织的优选法经验.这个科学方 


PDF 文件使用 11 p d f F ac 1 0 r y 11 试用版 本创建 www. f i nepr i nt . com, cn 







法已经 被千百万工人所直接掌握，他们应用这个工具，获得 
了数以百万计的优选法成果（这里指的是推广的成果数而不 
是发明创造的成果数).在国家不增加投资、设备、人力的条 
件下，为实现优质、高产、低消耗，怍出了十分喜人的贡献. 

这些年来：我们受到 r 极大的教育和启发，提高了对‘<群 
众性”和“实践性 t 的认识.我们体会到，一个方法如果只能被 
少数人所掌握，那末它所创造的劳动价值是很有限的.采用 
通俗易懂的平话形式，向广大群众介绍这个科学方法，是实现 
“群众性 & 和"实践性”的关键.但也有少数人产生了误解，认 
为平话没有什么理论.其实不然，从六十年代以来，我们就学 
习了几十种优选法，首先进行了理论上的推敲、剖析和研究， 
然后去粗取精，去伪存真，最后才挑选出几种理论上靠得住、 
又容易推广的和有效的方法.选定之后，我们又把这些方法 
深人浅出地逋俗化，这样才能被广大群众所掌握.实践 证明： 
在理论上深入固不易，而深人后浅出更不易.没有理论上的 
深入就没有方法上的浅出. 

我们认为，一个科学成果从理论上研究到把它普及到群 
众中去，其中差距是不小的.后者涉及到千百万人的实践，因 
此，.必须采取十分严谨的态度，必须经过反复实践. 

这本《优选学》提到了不少不适宜于普及的东西.例如， 
关于多因素问题的处理，可能不下于国外七十年代所提到的 
一箜方法，在理论上也可能更严密®.但是，我们还是应当以 
“抓主要矛盾”的思想为指导，在错综复杂的因素繁多的问题 
中，尽可能地抓住形成主要矛盾方面的因素，并且因素的个数 
也尽可能少. 

现在我们伟大的祖围正进人社会主义现代化建设的新的 
迈史时期，提高整个中毕民族的科学文化水平是实现新时期 
总任务所必要的，从这几年广大群众掌握和应用优选法过程 
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由促进了学习技术、提高业务的情况，反映了这一科学方法在 
群众中的传播与推广，可以为实现这一宏伟目标作出一些贡 
献.从运用优选、统筹听创造的大量份值中，也有力地显示出 
科学就是生产力这一马克思主义 原渾的 正确性_ 

在近年来的实践中，对如何把科学方法变成广大群众手 
中的工具，对如何洋为中用，对如何组织群众性的科学实验活 
动，以及如何发展适合于我国国情的科学技太等一系列问题， 
对我个人来说，总算与从前侷促于书斋之中，受囿于专题之 
下，不从国家需要和科学技术的相互关系考虑问題有所不同 
了.尽管自己在认识上提高了一大步，但是，这种认识还是粗 
浅的. 

这本书虽然推迟了十年才出版，但从学米水平来说 5 还 
未曾落后.为了把全党工作着重点转移到实现四个现代化上 
来，为庆祝国庆三十周年，谨以此书聊表心意，献给伟大的祖 
国，献给伟大的人民. 

最后，对参加本书校阅工作的龚昇、吴方、潘承烈、谢庭 
藩、徐伟宣、陈德泉、那吉生、裴定一、计蕾、李之杰、方伟武等 
同志，以及在工作实践中予以协助的推广优选法小分队的同 
志们，表示衷心的感谢，同时希望读者和参加过这一工作的 
同志继续提出宝贵意见，以便进一步修改补充. 


• m 
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第一部分单因素优选法 


It 


我们常假定/00是定义在区间（《， o 的单峰函数，但 
的表迖式是并不知道的，只有从试验中才能得出在某一 
点心的 数值外 nKO . 我们的问 题是： 用尽量少的试验次 
数来确定/0)的最大值的近似位置. 

最简单的方法是穷举法或称均分法，即在《个点 


Xi = a +— ~T ^ — 0 1 J -, n) 

做试验得^ = 其中最大的一个是 则我们知道最 

大值在区间 


O — d )， 


1 


(厶一 a ) 


中取，也就是做《次试验 M 以确定最大值在长度等于 



的区间 巾取， 如梁预先给了精密度 e ， 则试验次数〃必须适合 
于 

—-—(办 一 a ) S 


才能达到目的，即: 


ff + 1 ]> 


2( 占一存） 

e 


所需要做的试验次数的数量级是 i 才能达到目的. 

£ 

以后所讲的黄金分割法及分数法，所需要的次数的数 a 

1109191 … 
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级只要 log 1，大大改进了以上的结果，而抛物线法所需要的 

S 

次数的数量级 R 要 log log ^ s 又进一步改进了黄金分割法的 

6 

结果* 

就试验的次数來说，方法一个精似一个，但这几个方法都 
是有序贯性的，即必须利用上一次试验的 结果， 才能决定下 
一次做试验的地点，这样做可以节省试验的次数，节省人力、 
物力及试验费用消耗.但是，可能要拉长时间，如果‘件许 
可（例如有多套试验设备），我们可以用第三章所介绍的方 
法-批做几个试验的方法. 
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第一章黄金分割法和分数法 


述 要 


《优选法 平话》 一书介绍过 G .618 法. 0.618 是 

5 ^ 1 

w =--- -= 0 . 6180339887 —* 

的四位有效近似值，它适合于方程 

如 2 + 奴一^ 1 =» 0. 

黄金分割法也就是先在 




nam 


a -f* (b —■ a^)w 


处做一次试验，再在 n 的对称点 

x 2 = h —— (b • — = 


a 


O 一 




2 


做一次试验，比较试验结果 h « /(A) 及 y 产 f(x 2 ) 哪个大. 
如果 Kh) 大，就去掉 （a ， 在留下的范围 Oh o 中已有 
了一个试验点然后再用以上的求对称点的方法做下去. 
如果 Kh ) 大，则去掉在留下的范围中已有 
--试验点同样用求对称点的方法做下去，一直做到所需 
要的精密度为止. 


在介绍分数法之前 5 引进数列 

Z 7 。 = 1 ， Fi = 1 ， F„ = F ft _! -H F «_ 乂 /J > 2 )， 

即为 

1，1，2, 3, 5, 13, 21，34, 55, 89, 144,…. 

第〃级的分数法以 

F #+ i 

* 3 舞 
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代替黄金分割法的做了《次试验后，知道最 

大值在长度为 — 的区间中，因此，给了精密度 h 取曰使 

匕 +1 

F n+l > — > F *, 

€ 

则用第 〃级的 分数法做〃次试验就可以达到要求. 


S 1 . 有限点的问题——分数法 

不坊假定所考虑的区间就是[0,1〗，在仿个点& ^ 

0 <： A < : fj < : 3 <，■- < X ffl < 1 

处可以做试验 t 问题： 至少要做几次试验可以找到 lUi ) 

中最大的一个？或者我们至少要做几次试验， 
怎样安排，可以从最多的点数中找到取最大值的点来？ 

令 一武 是 ”次试 验可分辨的最多点数，也 
就是在个点中，可以做〃次试验找到一个^，使 
y _ = /CO (1 < i < m ) 为最大，而超过叱个点，我们就无 
法用《次试验找出取最大值的点了. 

首先，做一次试验0 « 1)， K 能分辨一个点，因此 A = 
1. 做二次试验 （《 «2夂只能在两个点中分辨，因此％ — 2. 

考虑一般的情况：在叱 —000 个点中已做了两次试 
验①:^②如下图： 

O O *■ ■ ■ O O 0*^*0 O O 1 ■ • * » o 

^—— . — ' —,—②—^—^ 

成 6 e 

即在《个点之后做一次试验①，再在*个点之 后做一 次试验 
②，还余 r 个点，则得 

0„ = + CO 

如果①的试验结果比③好，则最好值一定在 
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O ^ O 0***0 

%_« _^ QP —一 

a b 

中，连①在内，在其中共可做《 — 1次试验，可分辨的点数 
最多为 h : 因此 

a + 厶 + 1 < 少》-1. (2) 

又在前面。个点中最多只能再做 〃 一 2次试验，所以 

« < ^,- 2 * ⑶ 

同样，如果②点比①点好应当有 

办 + c + i < ® 0 - i , (4) 

综合⑴， 0),(3), ⑷，得 

1. (5) 

当且仅当（：0,（3)， （4) 均为等式时等号成立.即当安排方 
法符合条件 

a = c = ^>„- 2 C 对称性 )3 ( 6 ) 

從+占+1 1=5 办 + r + l »少逋-1 

时才能使可分辨的点数最多.这时 

巾， 一 u U 1. (7) 

改写成递归公式 

= 1, ^-2, 3>,+1 - ( ui )+ d + i ). ⑻ 
令 Fh+i = A + 1，则得 

™ 2, = 3, F n ™ ■P»-i+ Fn (« = 4, 5 ，…）， 

如果定义 F 0 = F x -1, 上面的递归公式可写为 

Fi = \^ F „ ^ Fh + Z ^- jO — 2, 3, - ■.). (9) 

结论：做《次 试验顶多可以分辨出 （ F n +1 _ l ) 个点中取最 
大直的那一点.安排方法是先 在第^ 点处做试验，再在其 
对称点（即第 F „ +1 - - F „- x 点）处做试验，如果左边的试 

验点好，去掉笫^点及以右各点，余下 （fn — i ) 个点，依法 
选行下去： 、 
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o o 


O X o - 


O X 0 O Hfe H » ii O 








如杲这 （匕 +1 -1) 个点在 （0, 1) 中是均 匀的， 即均分 （0 5 1) 
为等分，经过〃次试验，最后余下长度为■^一的区间， 

在（0, 1) 中取最大値的点一定在这区 间中.试验 安徘是 ，第 
一试验点在 FjF n+l 处，以后取它的对称点 3 这就是第 a 级分 

数法 琪一—片市 l 做 t « je •本 nr 
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M 






y 


便有 F ,- F a -^ F n -^ 

, 由归納法证明了公式 （1). 

这里的证明是假定了已知公式 Cl ), 再甩归纳法来证明 
的.如果不知道公式（1)，如何导出这公式？ 

令 

fix ') = F 5 + F^x + F 2 x 2 + ..■ + F n x a 十 *"， (2) 

(2) X ( 一 r )； — xf ( x ^) =■ — Fax —尸：/一 ... 一一 F n ^ x n 一 - * * , 

(2) X ( — r 1 )； —^ 2 /( x ) 一 F ^ x 2 —… .一 F rt 一 〆 —*•， 

三式相加，得 

(1 一 r — /)/(芏 ） w F 0 + O ' — F^x 

-H CF„ — F n -i — F fl *2) xn = 1 ^ 
m / GO ' 、 一 ” (3) 


令 

容易算出 


1 _ x — = (1 —— ^1^)(1 一 



O ) 


于是 

/ O ) = 




(1 一 似]:)(1 — U r 2 x ) 
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A s ( u^^y + a 2 






= f > 






S (次吋 + 


( 5 ) 


?? -■ - 0 


其中 


= lira /( 尤 〕(1 — « p iJr ) — Em 


ti^i 






w 


W 2 X Wi 一 W 2 


同样 


為 


似 1 


% 


(5) 式可写为 


W \ — w 2 




* 


1M ^ S 


« l 


ro ^% 


( 6 ) 


— 0 叫 一 w 2 

与 G ) 式比较系数，便得 

1 


F 




U/l — - U ^2 


( —如?七)， 


⑺ 


即为 （1) 式. 

这个等式告诉我们，虽然所有的 F fl 都是正整数，可是它 

5~ — 1 


们可以用一些无理数表示出来.由于 


当"是奇数时，则^等于 



2 



1 1 

"+ Vl"\ 

V 


\ 2 / 


»十1 


-< 1，因此， 
的整数部份， 


当《是偁数时，则等于 


1+\/^ 


M-H 




2 


的整数部份加 1. 


读者可参阅《从杨辉 H 角 谈起》 (青年数学小丛书， 中国青 
年出版社）. 

作为练习，请读者求出下面两个数贯的解析表 达式： 

1) 数贯由递归公式 

=1，+ m ) (^> 2 ) 

给出.当 / = 1 时，就是 


世 
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2) 数赏由递归公式 

r 0 == = 1, 5=5 c n -i c n - 2 (n > 2) 

给出. 

它们将分別在第二章和第三章中用到. 

I 

53. 黄金分割法 


* 


分数法虽好，但预先需要知道做多少次试验，这是很难 
的，因为预先给定精密度 G 吴差范围 ） e 是不容易的，一般说 
来，函数的变化不是均匀的，越靠近极大值，变化越平坦，一且 
试验结果不能分辨优劣时，试验就只好停止.所以，我们希望. 
能有一个预先不限定试验次数的好方法. 

现在我们来看看，当《趋向于无穷时，^的极限如何. 


因为 






1 


lim 


F 


73 


F 


- Km 


< 1，易见 

2 


«+1 



V 7 


#+i 


2 


fli+i 


W 7 

2 


tt +2 





li+a 


2 






2 


p 


我们记 






2 


进一步考虑差 


F 




议 a — 




F 


n+i 


F 


fi-M 
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«+2 


(一 iy 


F 


fi+t 


2 



5 


n+2 


F 


2 


»+l 


一（一 I ) 



^ \ «+l 


F 


1 + VT F 


•+i 


2 


2 




F 


l 


w 


F 


1 + \/ 5 \ f 

2 / f 


n 


( 


F 


F 


ji+-i 




F 


F 


A+if 


< 




•+i 


2 


F 


冪 


因此 


F 


F 


w 


ti+i 


< 


w 


F 


F 


而由 


F 


F 


w 


(— 1 ) 


n 


*-H 


I + W F ， 


n+l 


2 


易见 ； 当 《 为偶数时 ， A > « v 当 rt 为奇数时 ， i w 


F 


F 


^■ j-i 


于是 


> .^ > £< > 


Fj 

f 3 f 5 f 7 




心 < ^ < ^ < 

F 2 F . F t 




这说明 




F 


当 《 趋于无穷吋，在 w 的一左一右摆动，而且一 


n+J 

个比一个更接近％ 
用黄金分割法如下: 


10 
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I ， Xjt JTi 身 


再在 （h 6) 内取一点 
做试验得如果 

而则丢掉 
(〜6)，再在 Oi 5 中取 
一点〜 ■• •，不断傲下去，就 
这样通过来回调试，范围越 
缩越小，总可以找到 / O ) 的 
最大值（图1乂 

这方法取点是相当任意 
的，只要取在上次剩下的范 
围内就行了，我们的根本问 
题是：怎样取 A ， 〜，…， 
可以最快地接近客观上存在 
的最高点？也就是怎样的安 
排试验点的方法是最好的？ 
何谓 w 最好〃的来回调试 
法？先看一例，有一方法，第 
一试验点在2/10,第二试验 


点在1/10，对峰值在 （0, 
1/10) 中的单蜂函数，两次 
M ^ 试验便去掉了区间长度的 

4/ 5 (图 2), 但对于峰值在（2/10, I )的函数，只能去掉1/10, 


这就吃亏了 . 由于我们试验的目标函数是事先全然不知道 
的，因此我们说的"最好' 就不能对个別的单峰函数来说，而 
应对全体这类函数而言. 

这里不引入更多的数学定义.，只是根据上面所说的事实， 


在评判一个来回调试法的效率时，作如下的约定: 


如图3在< &) 做了两个试验以后进行比较时, 


• 12 * 
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I ， Xjt JTi 身 


再在 （h 6) 内取一点 
做试验得如果 

而则丢掉 
(〜6)，再在 Oi 5 中取 
一点〜 ■• •，不断傲下去，就 
这样通过来回调试，范围越 
缩越小，总可以找到 / O ) 的 
最大值（图1乂 

这方法取点是相当任意 
的，只要取在上次剩下的范 
围内就行了，我们的根本问 
题是：怎样取 A ， 〜，…， 
可以最快地接近客观上存在 
的最高点？也就是怎样的安 
排试验点的方法是最好的？ 
何谓 w 最好〃的来回调试 
法？先看一例，有一方法，第 
一试验点在2/10,第二试验 


点在1/10，对峰值在 （0, 
1/10) 中的单蜂函数，两次 
M ^ 试验便去掉了区间长度的 

4/ 5 (图 2), 但对于峰值在（2/10, I )的函数，只能去掉1/10, 


这就吃亏了 . 由于我们试验的目标函数是事先全然不知道 
的，因此我们说的"最好' 就不能对个別的单峰函数来说，而 
应对全体这类函数而言. 

这里不引入更多的数学定义.，只是根据上面所说的事实， 


在评判一个来回调试法的效率时，作如下的约定: 


如图3在< &) 做了两个试验以后进行比较时, 


• 12 * 
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图2 


图3 . 

总是遇到坏(吃亏)的情況，即只能去掉和 0中 
较短的那一段 • ' 

这样做了《次之后，余下一个包含有最好点的区间 ，其长 
度称为（该方法）《次试验的精密度. 

对黄金分割法，假定原区间校度为1，则做〃次试验后余 
下区间长度为 W 1 ， 于是 

因此，若预先给定精密度 S ， 这时所要做的试验次数《应适合 

^* -1 < 6 , 

即 

« > 4 . S log — + 1 * 

e 

也就是说 n 的数量级是 log 卜 

在所有的来回调试法中，以黄金分割法为“最好'我们有 
基本定理任何一个来回调试法 ，如 果不是黄金分割法， 
总可以找到正整数 V ，当试验次数《大于 JV 之后，其精度 
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总不如黄金分割法做同样次数 n 的精度 In ，即 

A n Iti = (n > N\ 


§ 5 . 黄金分割法的最优性 


在这一节里，我们将证明上节的基本定理. 

先建立一个不等式， 

对于任葸一个来回调试法，假定第一点和第二点分別在 
巧和~处做，不妨假定区间就是（0, 1), 且& <心，显然 

厶！ = 1 (图 4 乂 


A ^ t 

0 * 

m 4 

当&+& > 1时 ，心 > a , 不论第三点如 何取， 总有 

■ 

△3 > *1, 

故得不等式 

s C 1 ) 

当 A 十 h t 时 ， Al > i — h ， 

不论第三点如何取，总有 A ,> l - X ^ 

故 

Aj + A 3 > (1— A ) + (l — : 3 )二 2 —( r , + ^ 3 ) > 1=厶“ 

于是，不论第一、二、三点如何取 3 (1) 式恒成立. 

一 般地，我们有下面的不等式， 

△KH ^ 么 ^：+1 + Ak +3 (尺 ■’ )* ( 2 ) 

由 C 2)， 用归納法容 芴证明 

< F 入 + FaA tt+K +' ( 人 % » = 1 ， 2 ， . • .)• ( 3 ) 

这里 F 0 Fi ^ l , = Fn-i + 匕 -j ^ 2 ), 

记黄金分割法第 &次试 验后余下区间艮度为已知 

H ， 
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Ik 





K 


y 



(4) 




对于心，恒有等式 

^ (允=1，2，…）， （5) 

= Fn / fi+K 十 ^ n - L ^ w + K+l (方 ，人 ^ 】，2， * ’ .）. （6) 

下面，用反证法来证明基本定理. 

若定理不成立，则对于给定的来回调试法，必有一无穷数 

贯 < & < . ■ . < < …， 

使得 

A fljr < h s ( j ^=1，2，- ..)■ ⑺ 

对任意的 X ，必有〜，使 KCn ,< n i+[ < ‘ ■. 不妨就 


当作 « 0 -= A ： <«!<«,< "*. 

由不等式（3)，得 

^ ^» l -n () ^ J + Fwi(FwlA ni + F w 2 A n ，i) 




-H F 






■[ iL^h 






而白 (6)， 有 






付 j 


( 8 ) 



由于 < ’np △^K △ if ; < 



攀 15 * 
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从 （8) 式戚去 （9) 式，得 



(II 

U ^0 


F 




) 

/ 


F 




2 ^ n/+2 


do) 


因为 

F n — ( — 如 )…) 

5 


< 





由 Cio ) 即得 

Ax — /jc ^ . 浓一〜十 + 2 


w 


» 


•+I 


= uf s ^ n ^ 2 = (s = 1 > 2 ? ■. . )* 

令 f — O 0, 即得 

△K<U ( 尺， 1 ， 2 ,…）. 00 

由定理假定，给定的来回调试法不是黄金分割法，因 I 衍必 
有某一个6使 


A f < l f , 

由 (3) 和（6)， 

1 d A| € F + 

< Fs-ih + ^ \ y 

得出矛盾. 

证毕， 


S 6, 连分数的知识 

现在，我们换一个方法研究来回调试法，给出黄金分割法 
最优性的一个直接证明，并说明一痊更深刻的内容_ 

* 16 • 
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为此，需要用到连分数这一工具.如果我们约定只考虑 
对称取点的来回调试法，那末就与辗转相除法、连分数关系+ 
分密切，抽象地说，他们实际上是一回事、所以，我们先介绍 
连分数的知识., 

1〕先从两个正整数《，6 谈起.假定〜用 6除 a 
得商叫，余 r， 

a + r ， 0 ^ r <1 

若，今 0, 再用 r 除 6 得商，余即 

b 通 a x r -f- r iy Q ^ <C r; 

若 h 今 G， 再用 n 除〃，得 

r a 2 r x + r 】， 0 < < r ； 

若〜七0,则继续做下去，由此得 

办 >r>q>r 2 > …>« 

经有限步后，余数必为0 ,例如^ « 0,这个方法叫做辗转相 
除法.上面的过程可以 写为： 




这一繁分数称为连分数.为了节省篇幅，把它写为 


例1 





■> 
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般地有 

F 

■ 

F 


A+i 


Fr ^ 

F. 


十 1 


% 


1 个 


2) 辗转相除不一定限于有理数对任一无理数 A 也 

b 


可以写为 


of 在 d 0 + — ， a 0 ^ [ a ] , 巧 > 1， 


方括弧代表《的整数部分.再令 


， a t ^ [〜]， Oj > 1 # 

(*2 


一般地 


= a 


a 




[ 1 • 


因此，得到 a 的连分数表达式为 

, 1 1 1 

汉0十 ^ ― r * — 

u { + 这 】+ • • • + a n ^ x + • • ,• 

因为，任意有限项的连分数都是有理数，所以无理数的连分数 
必然是无穷的， 

经计算得 


汶0 + 


ao 


1 


fifj + a 2 


aaa^i + 斗 

a 2 ^i 4 - \ 


令 


<70 十 一 

A 




+ • • _ + # q 


y 


n 


称为 《 的第〃 个渐近分数, 


19 


» 丄 3 • 
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例 2 令 w 


7 ^ 




2 


它适合于 w 2 




1 = 0 , 


即 


w 


w 


W 


W 


■4* 


因此」 

F , 

? 

，就是 


F 2 



F , 

f 5 

f 7 

即得 




就是&的渐近分数，§ 3中的结果说明 




F ； <f ； <f：< 






> w > 


的一个小于一个的区阆套，又由 

jFj/ JF，;— I 1 


F21 


F 


ii-b 1 


Fa 5F 2i F 2i+l 


{[ 


w 




(-«/ 严 M 2 


5 尸 h * 尸 h+i 

1 

5 F ,, ^ F 


- [浓 - a — (—… — ( 一斯 ) 2 ^]} 

2 




^ + 2) 


U 〜一 (一 w )] 2 


Ai J 2i+1 

可知区间套一个套着一个. 




h % 连分数与来回调试法（一） 

考虑对称取点的来回调试法. 

在区间（0, 1) 中已做了试验的点为 

( i ) 若 ^ 上丁 ，0 < 0 i < 1* 

1 - T - c?i 

^ 19 * 
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图 5 

则《> 对称取点，第二个试验点在1 — «处，不论去掉 

那一段，余下区间的长度总是 
即 

Aj «« 1 5 A a oc « 1 

1+0/ 

不失一般性，不妨假定总是靠近左边的点好，也就是说，总是 
去捭右边一小段. 

因为 of 4- 0 x a =1，即 1 — a = B v a — 0^ 2? 

因而，从《开始做二次试验后，余下区间中的已试验点，即 
(1 — a ) 的位置恰好在区间长度的乂倍处，同样，若 

6i= rh-/ 

则用对称法取第三点，余下的区间长度为 

A 3 '= 0 iA 2 = 0 iOP , 

余下区间中的已试验点的位置在区间长度的氏倍处. 

■ ■ 

( iO 若 a -户为 >2的整数， 

p -b d 


困6 

则汉< 从 a 开始，对称取点做二次试验后，余下区间的 

2 

长度为 j 

A 1 1 1 P 一 1 + 沒 

△鮮 1 — a mm I — »- -~ = £ - - - 

p B p -b 3 

^ 20 ^ 
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其中已试验点(即 a ) 的位置在区间长度的 


⑴ 


OS ； _ Os) ^ - — -- 

— 1 1 + 0 

倍处，再做下去相当于在区间中从 A -^~ - 出发来回 

f — 1 + 行 

调试， 再做了 2 次巵，余下区间（ 0 , 1 一 （ p _ l »， 即 
由一段长为另一段为〃的二段所组成，所以再做一次，取 
点应在处了，余下区间的长度为心其中已试验点的位置 
恰好在区间长度的 0 倍处， 





0 a 

L _ [ ^ — " J 

t- iP-Utt 

图 7 

这样不断做下去，如 

ct = — 丄 - 不— ， 0 < < 1, 

A + + ， ■ . + -h 

则经 〜 + 次对称取点调试后(:不算《本身一 

次），余下区间中的已试验点的位置恰好在长度的0,倍处， 
再证明几个引理. 

引理 1若 。 … 




+ 1+6 


W 个 


0 < 石 < 1, 
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( 2 ) 



〆 


05 


P 


*1 


dF 


K 


这里 F 0 - 1, F 




证用数学归纳法, 


当 


1 时3 


TT ^， 


(2) 显然 成立. 


当対& 2时， 


a 


1 +9 


i + a 


(2) 也成立. 
当 m = 


时，令 A 


d 


则由归纳法假定，即得 


a 




F 


l -h 6 


F 


F 


1 + 沒 


F 


F 


6F 


F n + SF t 

敌 o ) 式成立. 

引理 2 不等式 

參 

1 一 > 


F 


fn+2 


5 u ? 


睿 


1) 


成立，当且仅当 


2等号成立，而 


I 3 + V 5 


5 u ^ 


10 


0.523 


证仍用归納法. 
当7» « 1时， 


F, 3 3w 


> 


5^ 


■ 


• 22 ^ 
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当 


2时， 


F, 5 W 


故当 




I , 2 时，引理 2 成立. 


F 


对一般的〃，若> %则由归纳法假设， 


F 


m+l 


F 


■jw -*1-2 


^ m+2 P«l+i 

r 


> 


( 忐 ) 


* 


5 


2 




, 则由 S 3, 


F 


F 


> 


即 ％±i < 丄，所以由归纳法假设， 


F 


# 


Fm+2 f F 


fW+-l 




F 


m+l 


F 


m 


> 




n / 



fi / 


z 


u / 


5 砂 2 


d 


故对■一般 《 引理 2 也成立， 

同样司证 

引理 3 不等式 

一 1_>丄 

P m +i 


>1) 


成立，当且仅当 


时等号成立 d 而 


— 0.745 




现在来证； 4 的基本定理. 
设区间 （h 1 ) 中已有试 验点久 

I .若 


* 23 


尤 ship 
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1 + - - * + 1 + $ 


, o<e m <\ 7 


iw 个 


a 


1 -h 6 , 






则再做扣次（共仍 +1 次）试验后，由引理 1 可知，余下区间 
长度为 


A 


相 +1 


C5 狀 … l-l K 

. ^ fn -2 


F 




+ e m F 






因为 0 < l < 1，故由引理3,可得 


* 


&抑+1 > 


F 你 + F 


F 




^ ~ w m ( m>lX 
2 sp 


H . 若 


a 


其中心 




1 -h " B # ^ 1 + dm 

n ■*'V —r ■ 

ra 个 

9 p > 2 ， 0 < 6 m+l < l w 


由于匕 < ★，酬 理 2 可得 


厶 


fKP +1 


p m + e M F 


> 


2 


F 




2 


F 


2F^+ F 


fii 




F 


m+2 


)W 


I 


m 会 1 乂 


24 


• « 


# 7 Bnhr 


i> 
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IH . 若 


Of 


p 0 


P 为 >2 的整数， 0<^< 1, 


分二种情形来讨论. 

0若3,则由 o )， 再做一次试验后， 


△2 


> 


2 

3 


2 


-识 ■ 


w 


p — I -\r 6 


2 


ii) 若 P 85=5 A (0, 1) 中已试验点在 ce 


ESi 


2十沒 


处> 第二 


个试验点应在1 — a 处，此时 ， a < — < 1 — 这与（0, 1) 

2 

中已试验点在1 一心 再做第二个试验点在《的惰形是一样 
的.而 


05 


6 


1 


2 +沒 2 +沒 


0 


归结为 I 的 情形. 


0 




L -« 


图 


因之，若 


ct 


p>2, o <e m < i. 


琴 ■ * 


+ 1 + P + 8 济 

在（％ i ) 中依 a 进行对称来回调试，由11，做 m + JL 次试验 


后，余下区间中已试验点的位置在 


P ^ 


， 2. 再由 III 


i )， 再做 f — 2次试验后，即共做 w + 1 + p —2 — m + /I —1 
次后，余下的区间长：度为 


• 25 
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再往下做就是 m ii ) 的情形了. 

现在我们回到基本定理的证明. 

若 a 是有理数，可由§8的论证中顺便推得基本定理成 

立. 

若《不是有理数，即《可以写成无穷连分数，此时，对称 
取点来回调试法可以一直做下去， 

如 a 七〜，那末《的展开式中一定不全是 I ,记这些不为 
1的数依次为九 ， Pa ph …•- 

I . 若这样的〜多于 < 个 (> 在下面定出），则当 ”充分 
大之后，綜合 I ， ii ， m ， 恒有 

[ 2 1 

2 

因子+= 0. 74 5 — _ < 1，至多是在最后一个试验点是属于 

Lti / 

1或 m ( ii ) 的情形吋才出现. ^ 

如果即 


S (八 _ 1) > 


log 2 




log 2 /W 


4,6 


■ 


时，则、> ln ， 只要取 c = 5就可以达到目的，即九■多于5 
个时，当 〃充分 大时， ^ 

A rt > 眇 ”- 1 — In. 

2. 若这样的 A 不超过 c 一 1个，则 


1 -4 - 1 -|- * * * -I- p c—l -f- 1 Hh 1 -4 - * * * 
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1 -\r m m m p { 1 + ■ pa + ti/ 


故做充分大"次试验以后，又变为黄金分割法了，此时 

4 ； 一 1 

A fl ^ I- 故卜 1 > 糾 0 一 1 in, 

\5 

这就证明了基本定理. 

附记衡量一个方法好坏时，还可以甩另一种精度，即 
以《次试验后余下区间中已试验点与较远的那个端点的距离 
^作为〃次试验的精密度.如〃次试验后余下的区间为（心 
办），其中已试验点为真正的最好点为〜， _： F 是 

1^0 一 ^ — a^b — x ), 

1 T 

« b 


图 9 

为了与前面所说的 A - 精密度区别，称这种精密度为 5 -精密 

度. 

分数法 的卜 精密度为心=+，黄金分割法 

^« L+I 

的心 精密度为 6 \ = 由于对称取点，总有心= △ 〜 .从 

上面证明即 得：在 > 精密度下，基本定理依然成立. 


§8•连分数与来回调试法（二） 

d 

上面证明了在不限定试验次数的情形下黄金分割法的最 
优性.现在再来看看用来回调试法限定做 2 V 次试验的情形. 
设 

_ 1 1 1 
a T + 了 + •. - + ' + & 

m 个 
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其中 m + p + lU 扣 

匕 +| <1作 




P + 


，整数0< 


a 



1 + 1+ … + l+(p — ]) + 


■ — v 

m + 1 个 

从 《 出发，用来回调试法 5 做 m 
的长度为 


次试验后，余下区间 


A 


tn+*l 


其中已试验点的位置在处，再做一次，因为0 < L < 


2 


,故有 


A 


+2 


(i 一 0 併 ） a 


— a 


nsHH 




p—i + e m 


+1 


(/? — 1 + UF 你 + F m+1 ^ 


其中已试验点的位置在 


1 


P 




A 


l »+2 


处. 

从/ 出发，用来回调试法，做 m + 2 次试验后，余下区间 
的长度为 


A 


m +2 


F 


m4-l 


沪—丨 + d m+l 


F 


m 


p- i-hO m 


+i 


^ + (f — 1 十沒料 i)^+i 


其中已试点的位置在 


P 


+ e 辑 


A ： 


m+l 


+1 


2 » 


暑 
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处，由于 

心 +0>— 1 + UL —&+■] 

口 （p — 2 十 6^+jF^ > 0 , 

故 


A 


* n +2 


> A； +J 


且它化 1 中间已试验点的位置是一样的， 故当 N > n>m + 1 
时，却有 

厶《 >厶“ 

或一样有 


也就是说，对于所有形式的 A 



£ t ^ + + ■ " * Hh Gtj Hr " 




其中〜， 々为正 整数 ， A + 心+ ■…+幻+ 1 = AT ， 

1 < AT — 1, 0 < 6< 1 # 


^ S = A Jl ―—， 0 C 0 < 1 

1 + 1 + ■… + 1 + 0 


出发，用来回调试法做 w 次试验所得的精度为最好. 

由此可知，只做 y 次试验的最好方法只要从&中去找就 
可以了，而从 s 出发，用来回调试法做 w 次试验后余下的区 
间长度为 


A fl 









图10 


* 29 * 
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其中已试点的位置在 处 # 
因而 






显然，当6? 


max (1 _ 6 ^ 0 ) A A r ™ max (】 一 8 、 

+ 0 F „- 2 


1 

— 6 

■^xV -" 1 

+ 

0 

_ 尸 N" - 1 " 

如 dF N ^； 


，若 0 <沒< 


2 


* 


2 


<ec 


* 


2 


时，〃 w 取最小值 


而此时, 


谷 N 


丄 

2 


F n - } + — F 


2 


N -2 


F 


JV+1 


# 


a 


iv - 1 个 


2 


i^=« 


F 


+ I F 


* 


w+i 


2V + 1 个 


这就是分数法. 


这样，我们在心精度下，证明了分数法 i 是 v 次试 

F 


A f +l 


验的最优来回调试法，这与 s 1中的结论是一致的 


显然 


衣 N 




F 


即 


N +1 


) < 占 JV («0, 


F 




W +1 


由引理3,可推出 


• 30 * 
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从而有 


即在心精度下，分数法比黄金分割法好.如用黄金分割法再 

多做一次试验就比分数法好了. 

在 A - 精度下情况就不一样了.由于如是0在（0, I ) 

中单调递减函数.当 0—1 吋，单调递减于 

1 — 1 
Fjj-i + Fn-7 

因泠只限于做 2 V 次试验，故0不能等于1，也就是没有一个方 

法能达到精度即不存在一个最好的方法. 

Fn 

分数法相当于6中0 = 1 -的情形，而黄金分割法 

尸訊 2 



相当于 B = ^ 的情形，由于 w ， >丄，可知 

2 

即 

此式亦可从引理3直接推出，这说明在 A - 精度下，分数法又 
不如黄金分割法了. 

到此，分数法与黄金分割法的关系就比较清楚了，并 a 解 
決了 a 是有理数的情形. 


S 9 .对分法 

上面介绍 的是： 求函数 /0 O 最优値的一些方法和它们 

- 31 » 


PDF 文件使用 11 p df F ac1 0 r y 11 试用版 本创建 www. f i nepr i nt. com, cn 



















的理论根据.在实践中也常遇到这样一类特殊的问题, 例如： 
在保证质量的前提下，如何降低某种原材料的消耗？如何降 
低成本？如何尽快地找到电路、通讯线路、下水道、输油管以 
及机器、仪器的故障？砂轮、转子的静平衝等.这些问题的特 
殊性在-丁_:每次拭验的结果只有二种性态，因此可以化为另 
一类的数学模型. 

对这类问题，我们建议用以下的对分法，而不用 Q .618 

法. 

例如，有一条10公里长的输电线路出了故障，在线路一 
端/处 有电，在另一端 B 处没有电，要迅速查出断头所在位 
置，我们先在线路中间的 C 点检查，如果有电，说明断头不在 
AC 而是在—段，又在中点0检査 s 如杲无电，说明 
断头在 C £) 部分，再在 CZ ) 中间 E 处检査，如此类推，很決就 
能找出断头位賈. 

i _ f ? ? f 

有电 有无无 无电 

对分法最简单，但在实践中应用非常普遍.用对分法每 
做一次试验，区间长度缩短一半，若原区间长为~做了 〃次 
试验后，余下区间长度为 



易对分法所需试验次数的数量级也是 

log — , 

8 

其中 s 是预先给定的精度/ 

对一些不易求解的方程或方程组> 用对分法能迅速求出 
近似解.可参阅本书附录二， 
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第二章抛物线法 


述 要 


不管是 （1.618 法还是分数法，都只是比较两个试验结果 


的好坏，而不考虑已做试验结 
果的数值如何.能不能 利用& 
得的试验数据，做到对技末精 
益求精呢？例如在试得三个数 
据后，我们可以过这三点作一 
抛物线，以这抛物线的达到顶 
点处作下次试验的拫据.确 
切地说是 这样： 在三个点&< 
a < ^，分别试验得数据 yi ， 
yiy 我们用插人公式 



O — x^)(x — 

— ^ 0(^2 — 文 1 ) 


y yi 


— _ 丈 3 ) 






O — 4)0 — xj ) 

— ； oo 3 — a ) 


得到一个二次函数，它的图像就是一条抛物线.这个函数在点 

龙4 =丄 一 ^0 + yz(xj — : j ) + 

2 ^iC^2 — 尤 3 ) + — 尤 1) + >sC^i _ 太 2) 

处取最大值，因此，我们下一次的选点取 r = r 4 , 得试验数据 
y ^ s 然后利用 d ， 〜， h 三点的数据乃， h 又可得巴一条抛 


物线，如此继续做下去. 


* 33 ^ 
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当然，为了更快接近最大值，可在 y 、， 仏〜 h 中取最大 
的，所对应的^作为新的并取其左边最近的^作为新的 
右边最近的; r , 作为新的再同样做下去，不断获得新的 

y 值. 

如果发生 A p r 2 的情况，我们的方法还必须修改 （5 1). 

粗略地说，穷举 法作〃 次试验的同样效果，黄金分割法只 
要数量级〃次就可以达到，本章提出的抛物线法，效果更 
好些，只要数量级 logh 卟次 ，原因就在于黄金分割法没有较 
多地利用函数的性质，做了两次试验，比一比大小，就粑它丟捭 
了，抛物线法则对试验结果注意到数量方面的分析.当然，要 
睜证抛物线法收敛速度，对函数也作了一些适当的限制. 

抛物线法常常用在 0 ^ 18 法或分数法取得一些数据的情 
况，这时能收到更好的效果.此外，我们还建议大家做完了 
0 * 618 法或分数法的试验后，用最后三个数据按抛物线法求出 
〜 并计箅这个拋物线在点 r = A 处的数值 5 预先估计一下 
在点 A 处的试验结果，然后将这个数值与巳经试得的(在 yn 
y 3? y , 中）最佳值作比较，以此作沟是否在点^处再作一试验 
的依据. 


S 1. 第一种方法 


设在三点 

C 

蚍分 別测得> = fOO 的函数值 

/C^i) = K 上 2) 二力， Ka) == y 3 , 

用 Lagrange 插入法作拋物线 


g ( flf ) = y x 




O — 々）0 — ^ l ) 

一 aXa — ^l) 


► 54 • 
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— 一 

■ ■■■■■ ~~ 1 1 l 

c ^ — — 


这个函数的极值由 


g \ x ) = yi 


c 


Jf 


2r — >y 2 — x 3 
— aX 吣一 a) 


> r 2 


2x 


丈 3 


A:, 


0 2 — Xi){x 2 — ， r 】） 


； Vj 


2x 


辜 1 


^2 


0 3 — X t )(x^ — x 2 ) 


0 


给砧，也就是 


X 


尤 4 


E=i> 


2 


尤 2 


=丄 ViCd — + y ； C ^ — jcD + y s Cd — ^0 

2 y! ( 巧一 ^0 + yiC^j — a) + yhi — 尤 2) 

(^ jdjCy? —— yi) + C^i 一 ^0( 少 3 —— yi) 

C^ 3 —— ^)(^2 Jl) 一 ^ 2 )(> P i — ^) 

丄 Ch — ^) Ky 2 — + Ci — ^ i ) Ky ^ — yQ 

2 C^s — x i)(y2 一 yt) + (^| —* 尤 2)(ya 一 ^) 


现在我们假定 


>2 = K ^ l ) > 


/ C ^) j 是向上凸的，也就是 

(A — 尤 .2)/(6) + (:2•— x l)K^) 


CD 


? 


A 


^1 


即 

(h — + (文 I _ 0乃 + (h — x ^) y s <0, (2) 

这时候 /XO < 0,所以 A 正是 gM 的最大值点. 

于是，立刻建议出第一个 方法： 如果 


yi > yi > ya 3 

那末，我们虫公式 （1) 得到在 ； t ^处做试验得函数 
y = /(^) 的值 y 4 . 若 > > 乃，取 h 作为新的 h ， 并取其 
左边最近的^作力新的右边最近的 A 作为新的若 
>4 < 仍取 h 作为新的再同样做下去. 

这样 ，我们 就面临着两个 问题： 

1) 这个方法做得下去吗？如果 A = A ， 怎么办? 

2) 这个方法能否导致/00的最大值？ 


• 35 * 


PDF 文件使用 11 p df Fac 1 0 r y 11 试用版 本创建 www. f i nepr i nt . com, cn 




















先回答问题2)，由（1〕式及（ 2 )式得到 

^4 = I (A + 文 2) 


— 1 — jQOj — X ^( y 2 — V 3 ) 

2 (% — h )( y 2 — yi ) + (^2 — -0(乃 一 yO 

< 丄（方 3 + A) ;> 

2 

和 

A =冬 (>■ + ; r 2 ) 

2 

+ 丄_ 0 3 —々）0、一 fa yi ) _ 

2 (^a — ^(yT^ yi) + ( 义 2 — ^i)( ； V2 — y 3 ) 

> 丄 （A + 文 2). 

2 

假定 

x t X 2 <i X 3y 

如果 M > 仏则函数 / oo 的最大值所在的区间由 o ls ; r 3 ) 缩 
成为（^2, ^3); 如果< >2,则缩成力 Gl ， O , 而 

X A ^ — (^3 + 尤3乂 
2 

同样处理 *4 < A 的情况. 

总起来一 句话： 每做一次，区间缩短一次，从而得到一系 
列的前一个包着后一个的区间 

W °，4 0 )， 。-=1，2,…）， 

因而其极限存在，记怍 

W '4 0 )) 

(请注意汁> 是不降的， ㈤ > 是不升的），不难证明4% 之 
一就是我们要找的点此处，我们不作证明的原因有 二:① 
留给读者作练习，@下一节所谈的第二方法的证明比这个更 
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深刻. 


现在冋到问题1)，当 = A 时，可以看看 A — :^和 
心一^哪个大，如果 A — A 较大，訢的点 h 可以按下法取. 
⑴取 

x \^ x % — 8. 

这儿 e > 0是要求的精密度，然后再用抛物线法做下去.如 
果 / GO 与 Ka — e ) 不易分辨，则说明 A 可能已到顶了，换 
取 W ^ x 2 ~h B t 如果此时迅是 /(^ a ) > /(rj 4 - e ) 3 则下一 
试验点就取在 h + s 2 的地方 等等. 

这样取点法优点是快，缺点是 e 难以预估，试验结果可能 
难以分辨+ 

( H ) 取 

」 _ A + 太 2 

" 4 一 - ， 

(若 h — X 2 较大，则取这顼法速度慢，假如适 

逢我们的目标函数就是顶点在 A 的拋物线，那就吃亏了. 

为了克服 （0 和 （ ii ) 出现的缺点 ，可以 采取 
( iii ) 在 （ i ) 中令 

S = ^(^3 — ^l) 2 

其中々为定常数，要求其保证无论进行到哪一步都有 

々(>3 — u) a < — (^a — 

2 

可令々为原始区间栓度二倍的倒数. 

然后取 4 — ^ + e ， 再按抛物线法做下去.这时显然有 

这样，既照顾到分辨问题，又照顾到收敛 

2 

速度同 A ¥ ar 3 的情况一致. 
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s 2 ■误差佔计 


假定函数在 r 处取最 大値， 抛物线法就是建 

议用上节公式 （1) 算出的6去近似地代替％，本节将紂误 
差1:4 — 作出估计， 

显然，如果 々不在 区间端点 5 则 

fOfl ) = 0. (1) 

由上节公式 （1) 可知 


方 4 — A -0 



(. r 3 4- x z — 2 xq ) 



记其好子、分母分别为/ V ，£)，并令 


= 

则 


尤 D + f _■ 


(i ^ 1 3 2 , 3 , 4 ) + 


N 


爹3 — It 


— 欠 2 






D 


由 T 


2 


Si 


r ： 


?2 §3 — lz 


A — y : _ y f 3, — y 2 

ll — ^3 — ?2 


一爸 2 


所以 


y ■- — ^ = }(x a H- I,) — /(.V 0 + I-) 

=(L — DJWo + 左 j .+ d — 


D 


2 


1 




(/ (^o + i — 0 + SiO 




— 2 


f (、 XQ + §2(1 — j 0 + &/)) 心 


ro 


/ 〃（方 0 + §2(1 —— /) + 


专 3 (i — 


4 


3S 


i-H-h-M&r-m I -^x- h* 
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I 


如果 


I 广 o ) 丨 > 吻， 


( 2 ) 


则有 

I 乃 I ^ m 2m (3) 

事实上，我们在上节就已假定曲线 _>• == fO ) 是向上凸的，因 
此（ 2 )含侖广 G ) € — m 25 于是从 Z > 的表达式得到 (3). 

再考虑分子，类似 _ f 上述，有 

n 口 —(Mi + ^)(( fX x o + 芒 2(1 — 0 + 

J 0 J 0 

+ &0 — ti)t)tdfidi + f f(x a -{- 6(1 — z) 

Jo 

+ 以)丸 
但因 fUX 所以 

fO 。+ 占 2 (1 — 0 + = f ( f(^D +6(1—0 

* U Jo 

+ ev ) — 厂 ( a ) V ， = H 厂(方 0 + — 

+ ^^)(^(1 — /) + § 3 ， V ㈣ ， 

及 i 1 

(\ ( Sz( 1 — f) + ^lOdptdi = 丄 （h + f 3 )， 

J0J0 2 … 

所以 

^ = — (h + 5) f (广 （A + 石 （1 一 0 I 

. 0 jo 

+ — 1^)0 — fX x ^))tdfidt 

⑴ 
⑴ 

. * 39 * 


+ 



0 


(广 （h + f 2 (I — f)/i + ^// i ) 


—厂 ( h ) X §2 (l — /) + ^ t ) dfldt . 
这样，如果 irwi 有界：即若 


irooi < 々， 


就有 
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I J Vj ^ j^3 -f- ^ 2 | — /) + 友冲 + 占 3(1 

JoJo f 

— fz)t\tdfidt 4 - | 0 | 0 ^ —0# 

+ 茬 wj ! ?^(1 一 /) + ^ \dudt 

十 fil (丄 ^ + 丄 ^ + 丄^) 

+ 《寻十 乎 X — 

+ ii + jg^l + l^il + \ Md \ 

3 2 - 6 6 y 

因此，我们得到以下 结论： 如果％ >0, 

roo< —叫， irw I < 々， （ 6) 

則有 

UM <^ max (! f 7 IS \ SA \ I 沾 h l^i t , I 沾 I )，⑺ 

Zftt 2 


S 3. 第二种方法 


根据上节最后一个不等式，我们可以求得另一方法，先用 
其他方法找到 


尤 1 < A < 怎 3 ，太 3 


lk 


及 

y2 > y 2 > yi . 

这儿 0 < l , 因此 

1^1 (; = 1，2,3)， 

3 弋 . 

现在用以下的方法来确定试验点，当有了心^， r …和 
h 之后，由 

40 


⑴ 


* . 
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[C-^n—i'~C-^-^ ~" 1 

+ C^n _ z"~ 1 iii]/[ C-^a - ' i - ~" i 

+ ( x n — ( r*-i — ^- i ) y n ? 

定义6+1，再测出 y «+ ij 不管 yn-ii ym ， 和 j ^+ i 中谁大谁 
小，我们把同一方法用在 〜 +1 上，下面证明这方法 
—定收敛，不仅如此:我们还要证明 


这里 q q = r 3 = 1， 

^ 〜 -5 O = 4, 5 …）. 

由 S 2 不等式 （7) 得到 


( 2 ) 

(3) 


Ifni < 弄 max (j Uh 丨 ， I Un ! ， 

2 tfT 2 

(4) 

因为 ^ 是一个不降的整数列，所以 

* 

- min (g- s + “- 3 , f•-! + ■+" —s ， 

2c n -i) — c n -i + c^-j (5) 

由 （ i ) 知道不等式 （ 2 > 当力 _ 1， 2 , 3 时成立，为了完成 
(2) 的证明，只要用数学归纳法就可以了^事实上，假设不等 
式 （2) 对 n » 1，2,…，々一 1都成立，那么由00, （5) 两 
式得到 


< 


3々 / 2 m ； 


2 m 2 ^ 




m 


in 公 Jfe-a* 


2m ： 


f 々一， 十【卜、 


即 


1^1 < 


2m, 


代 
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M . G 的表达式 


现在我们研究适合于 






fi-ly ^0 = = ^2 


的 G 的表达式.先算出几 项来: 


3 * J 


1，2, 2, 3, 4, 5, 7, 9, 12, H ， 2!，28, 37' 49 

斌 r * 的方法同第 - ^章§ 2 —样，令 








由 


cc 


(1 


JT 


一 /) S 


BC 


+ I + X 2 + ^ 5 ^ 




w »+3 


^ = 0 


ff»= 0 


40 




—— S ^+i: m+3 — ^ - ，十 j 


c^x 


f#= a 


/ —j 1 - #n 

o 


X 


得到 


/OO 




X 


1 一 jr 2 — ； r 3 


令 


i — y — / = (1 — w^X\ — 浓， )（1 — 如 〆) 

即％是方程式 

〆 一少一 1 七 0 

的根，这方程有一个大于 L 的实根 
对值小于1的共轭虛根. 


3 


A = L 324 •… 及两个绝 


用分项分数法， 






1 — U^iX 




為 


似 2 尤 1 — ti^x 


⑴ 


这儿 


• 42 * 
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A i = 


lim( 1 — 

h - 1 




(1 + x ) ti ^\ ± 
2x 4 - 3- r 2 I %:= ~i 


(1 + u f i)wi 



. 


由 （ l ) 式得出 


因为 

所以 


c n = A^i + Aztp^ -h A^w ^ 4 

limwj ^ lim = 0 

« -+K fT+CC 


于是，我们有以下的收敛结果 
即《又升高了一个台阶. 


附记在考虑抛物线法的收斂速度吋，如果把条件 
iro ) 丨 < 々换作 

irco - rv )! 

或者说 /" OOeLipU 自然是可以的.或许有的读者会提出 
roo < 一 不易验证，那末把它换作初始三点 
上 f ( x ) 的二阶差分 

—^…_ 之 L 二^) > 叫 > 0 (2) 

^3 一 ^1 ^1 — X 2 x 3 _ x 2 / 

也是可以的，这就表明，对于给定的/，开始用其他方法安排， 
一 旦条件 （2) 适合时，再用抛物线法，收敛速度就更快了，例 


如，此时若〜一: t ， < 


2 m 2 


‘ o <《< I ，就易于证明 


& 一❽1 


< 


2 m 2 

3々 






G O Q = 1， 






r =. 


4,5 ? 


* ■> ， 


乂 



# 
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S 5 .第三种方法 


从三点 Xi < x 7 < X 3 出发，有了 Jfl , ■ MJT * 之后要去 
找^ 在具体应用拋物线法时，除前面讲的两个办法之外， 
我们再建议用如下的算法： 

1) 先从三点 X x < ^出发 3 并假定 h <竹， y 3 < yi ， 

即中间髙，两头低，定义 

A ^ x 3> ^ h ， a^C c^C 

2) 有了 a ， 〜 ■ "，^ > 之后，假设 yi ， yi ，'' w y ^ t 中最大 

的一个在^ 处取，左边较大的一个在 X 成 a 节 处取，右边 

较大的一个在 ^ = K 处取，即 

〜 /(〜）>/(〜)， f ( c a ) > 

3) 令 y is > 2 , * ■，: y fl 中最大的三个数为 

Vi ** KD > vi 古 f(D > f(S^) 

由 （ 2 ) 可知 I — h 是 a ，< 之一，6或 是〜， h 之男一 

个，或在〜之左 或在心 之右.由这三个值作 

△?3 =专3 — 

Ah - — (^ + 6)， 

2 


△彻 取— V 2 , 

△ I2 = —吞 1， 

Ah -丄 （ mo , 

2 


A ?]! = ?]2 — ? J 1# 

再算出（参看 s 1，⑴) 




A ^ 3 A ^ Arj 2 ^ ^ g z Ag 3 Aj 73 

厶取一 厶1 厶也 


⑴ 
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4 ) 〜的定义如下 

f 若 h < d +1 < An ， 

X ^ _ 6 ，若 r :+i = G ，— a) b„ — c n - 
W+i + £ ，若 b n — c n > Cn — a n . 

A+i a { i 

了 ( a rt + ( n )， 若 € d n . 

+ 0„ + 心)，若 rD 

试验出 = /(A + 1) 的值，归纳法继续进行，这里 e 是个 

适当小的数，但应使 rf , + e 或 4 +1 — e 与 4 +1 的试验洁杲 
有所区别. 

附记0由以上的理论研究可知，到充分多次后，实质 
上总是 ^ : :4 a * 

2 )公式 （1) 的几何意义如 下：令 


即得 


X = ^^ Ag 3 A / j 2 — A ^ 2 Ak ^ A 7 j ^ 

^^3^372 — A^ 2 A^ 3 


X — M 


3 


△袍 功 y 


AT — 
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左边表示通过 d ( f 3 + 套 2 ), 0)且 平行于 小〕，（匕，… 

的联线的直线，右边表示通过 ( i (^ + o ) 且平行亍 

( S 2, ％)，（爸3,切)的联线的直线，这两条直线的交点的横坐标 
就是 -^*+1* 

这就是从抛物线上三点找出顶点的 方法. 


* 4 € • 
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第三章分批试验及其他 


述 要 

在生产和科学实验中，往往有这样的 情况： 做完一个试 
验要较长的时间才能看到结果，若等上次试验结果出来后，才 
进行下次试验，则要花费较多的吋间*如果试验还要进行多 
次，那末消耗的时间将更多.为加速试验的进行，常常采用一 
批同吋做几个试验的方法. 

例如，一批做四个试验，共做两批，我扪可以用如下的三 
种方法 安排： 

1) 00*00^00*00*00 

2 ) 00 0 ^ 0 ^ 000 ^ 0*000 

3) 0000**0000 朱 *0000 

在打* 的四点做第一批试验，然后在较好的一点附近四个打 
O 的点处做第二批试验.我们看到第三种方法共区別 16 个 
点，可以区别韵点数最多，但由于两对试验点的距离很近，试 
验的结杲可能不易辨別好坏 s 甚至会出现这样的 情况： 由于 
观察的误差，好的一段反而被去掉了.第一种方法少区别两 
个点，但第一批试验的点分布均匀，可以避免上述情況，昕以 
人们常常采用第一种方法. 

对于每枇做偶数个试验的安排是类似的.如果每批试验 
的个数是奇数，则安排便有所不同了，例如 ，每枇 做三个试验, 
做两批，也有三神方法： 

I ) ooo * * ooo 本 
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2 ) 000 本 **000 

3) OO 本 O 木 OO * O 

它们都可以区别九个点，但每种方法都有其特点，与每批做四 
个试验的情況类似，第一种方法的第四点、第五点试验结果可 
能不易分辨，而第九点则没有充分发挥作用，如果这么靠近 ® 
可以分辨，那末还是采用第二#方法好，因为在第一批试验 
中，如果第五点的结果好，它便是最好点，省去一批试验.第 
三种方法是为了避免出现试验点太靠近而考虑的，不过，如 
果要考虑做三批以上的试验，则第一种方法又是常采用的, 
从上述二例说明，在一批做好几个试验的时候，必须根据 
实际情况，选择适当的安排方法. 

S 1. 试验方案优劣的衡量标准 

在衡量分批试验方案优劣时，往往“一般地 35 采用留下区 
间的长度或最好试验点距留下区间端点的最大距离作为误差 
(或精度乂使之最小.在前面的傻定下，我们在 to , 1] 区间 
中点附近相距极小的两点做 试验： 


图13 

便可去掉区间的1/2—匕于是，也就不可避免地引人“在同一 
点作两个试验可以去掉区间一半”等这样纯数学的假设 • 

我们从实际应用出发，要求：最后一批试验点和上批留 
下的较好点等分上批试验后所留下的区间_这样使得全部试 
验结束后，听找到的最好试验点，距留下区间端点距离相等. 
在上述前提下，我们用最好试验点与留下区间端点的距离作 

_ 4S _ 
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为误差，并使之最小. 


§2- —组特殊的方程组的解法 


给出方程组 


〆 


£=sa 


议 0 == 2 1 « 1 + fiia 2 y 

(h 一 十 f«Wa ， 


⑴ 


a l-l 又冲 + ㈣ 窃 PH9 




^ 句 +1 


其中 - ，，， V ， • ■ * ，灼为已知数，〜，*•.， 

^ + i 为待求数. 

此方程组固然可以用最后一方程代入上一方程，解出 

然后又代人前一方程 i …5最后解出 
此方法较烦琐，我们可采用如下巧妙的 方法： 

令 

^ - =* (2) 


a 


众 nHKl 


则⑴化为 




X 


«1 


^1 
於 1 




1 2 + -^ 
^2 


友 2 


^ 


(3) 




l t + 

vi 


n 


由 (3) 3可将…〜…^ W 表为连分数: 


* 


49 
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Pi 


令 


IV-i = ll 


1’ 卜 2 = 又 / 一; 




I 






A 




/ 


«/ 




«/ 


^2 = ^3 


户 3 丄只 4 


足 4 Jt 


妁一 


⑴ 


5 




1 +苎+丑+ 
-^3 又 *4 




/fl 

又2 


T , 


Xl nr / i / 
+ 只 / 一 ] 

h ~ f - m 


^i 


m 


Fq — F 


则由⑷得 


?rt 


m ~ ^ y 2 y … ■，/) 




C 5) 


由 （2) 和 （6)， 


= 

= 




f 2 

^-1 

■— W - r ™ 4 

L ■ ■ 



斗 2 

f 3 

=._? 


尸 2 


F f , c 

^2 - 



F 卜 1 

Vl 一 

Fi 



^0 =*= 



CO 




a 


k 


^1 F t 


F 


V ^* Vi m ' * JJ 




F 




出 F f 


F 


/— 


50 


■ 


rHM -frH •中 V % HMfWmb -fj.. 
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I 


Fi 


i±l 


F 




= + 即有 


^0 

1 — 


m ■ 

■■ 

■ ■ 

■ 

a l 

Fi — 

h 


■■ - 

■ ■ 

a k ^ 

』 F 卜 
F I+l 


⑺ 

a J-i : 

f 2 

^^5 r - - ■ - - - — 

^ F 2 ^ai 

■ 


A 


F 




f+i 


F t+l 


§3. 每批作奇数个试验，如何安排 

试验 范围〜 =1，共做/批试验第一批七个试验，第二 
批 b 个试验，■…，第 J 批幻个试验 5 其中••_，&均为 
奇数.试验妄棑方案 如下： 

把別分为长度为 A ， a 2 Ca t > ai ) 互相间隔的（七+ 1) 个 
小区间： - 



m h 

试验后 留下包含取较大值试验点的 区间: 



图 15 


* 31 * 
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第二批 心个试 验全部安排在长度为 A 的区间 里，并把 A 分 
成长度为 b 互相间隔的（七+ 1) 个小 区间： 

…^2^3心2或 <?2 ^ ㈣ 2,， •这 钟1 

试验后又留下包含取较大值试验点的 区间： 


㈣ 3或 ¥2 

把下批试验全部妄排荏较长区间里，并把其区间分为长短相 


间的 Ob + 1) 个小区间，等等.根据 s i 的要求 s 第〗批的心 
个试验点要将分为与^相等的(:七+ 1) 个小区间 5 最后 


留下 



m 

此处的^即为我们所定义的误差. 

把上述的{^}的关系写为数学表迖式，就是纟2方程 
(0 的形式，其中 






2 

气，土」 

― 2" 


U 2, 


* 4 > » 


，0, 


⑴ 


0« _ 1，2，… W )- 


根据 S 2中 （7) 式可以算出而误差4 
面的递推公式 确定： 


F 


由下 


i-bt 


F 


心〜1， 


F 




fit-M 


2 


( F w + 1 ，2，…，/乂 （2) 


下面举几个特例： 

例1只作一批々个试验. 
由（ 2 )式得 
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Fq ^ F i = 1 ^ 1^ = 是 +1_ 


再由 S 2 中 （7) 式知 


A = 1 j 


a 


f 2 


+ 


即々个试验点把区间 A 分为（彳+ 1) 等分，试验误差 


这 2 


辱 


例2共作/ 批 ，每批作一个试验 
由（ 2 )式得 


F n = F 


F 


irt4+l 


= u 

尸 m + f m 


(m = 1 ， 2， ...，/)• 


这就是通常的分数法.试验误差为^ 


F 


* 


i+i 


例 3 每批作 （2 i — 1) 个试验，共作 々批. 
这时 （0 式变为 


Fn= F 


F m +i = + m = 1 5 2 > - *. ，々 )_ 

m 笫一草 s 2 的方法：可得 


F 


nHH 




2 


+ 



2 


— 3 


t 


•V 


h 4 


\Jt^t - {- 4 ) 

2 

— V 才（, + 4) 
2 


m 




且有 


a 


Ia 




F 

F 


Jt+i 


k-i 


^+i 


往 k j 沒 jfc+i 


F 


m 




*53 


(3) 
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试验误差为 




F 


* 




4. 每批 作偶数个试; 验 ，如何安排 


设笫 一批幻 个 试验， 第二批 h 个试验，一，第 〖批 6 个 
试验 ，七， 々2,… •，七 皆为偶数，把试验 区间〜 (=1〕分为长 
度是 〜 +1 (〜 > 句 +1 ) 互相间隔的 〔I + 1) 个小区间： 




a 


x—x 

人-入 




:X -X- 
人 J 

U.^i 


Tv 




^1 


17 


第二批试验全部安排在长为^的小区间，并把^分为长度 
a 2 ， a w (^> ^ +1 )互相间隔的仏+ 1) 个小区间，等等，根 
据 S 1的要求，第/批的七个试验要将区间分为与 
相 等的如 + 1个小区间： 



图13 

此处力 +1 即为我们所定义的误差， 

把上述的 {〜} 的关系写成数学表达式 




2 


A + f 一 


^2 1 \ 々2 

a 3 == I -h 1 } a 2 ~h — ^； +1 . 3 


⑴ 


1 


2 


〜廿 ^- a i+i 3 
2 


这 i ™ ^i+u 


54 


• 、砷畢 
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或 


这一 1 + a I+l 





( 2 ) 


所以 


1 + ai+-i = 外 + 这 h«i 


k 


2 



z 



(A + A+lX 




2 


k i 



. 


+ 


(3) 


再由 


a l+i = ❿ + A+i 



)( 警 



(〜 + A + i )， 


可以算出〜（他=1，2, …， O . 


a 


m 


々 pa+l 一 u hjn^l 


2 



+ 


2 



h 


l 一 


(4〕 


例每枇作以个试验。共做/批. 
由 0) 及 （4) 得 


i 


^ i+i 


2(j( + 0’ 一 


a 


2(主 + 1 ) 


/—m 


m 


2(/ + iy 


學 


55 


■ 響 
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§5. —般情形，如何安排. 


和5 3, § 4方法一样，若第 — 1〕批后留下区间为 

I ^ 癱 ― 1 t a P I 


图 19 

CD aO , 则第 m 批试验点全部安排在较怯区间〜 - i 内， 
当^•是奇数时，4 =‘；当 b 是偶数吋， a p = a m+im 
是 {相+1 > 1+2, … 中出现的第 一 个奇数 ， 

具体地说■•，，&中的奇数依次 为狄， 々 y ■， 
i it , 则各区间长度的数学表达 式为： 


1 — f?{] :== 

A + 1 

汉 j Hh 1 

V 


l2 」 


l2 J 




« 55 * 
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同 S 4 —样，可以把偶数次的试验合并写成 






fe 






一 贫、 9 









2 


k “ 


+2 


h 

L 2 


J 


2 


(A + a +1 ) — a M ， 


沒 i+u 


⑴ 


进一步可简化成（记 = O 

f 1 tfp = -h 

^ i L = ^ i t + 

^ . O ) 

5=1 ■+ 

、 ai = a t +u 

在这里， 
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I ' 


2 




i +^ 


2 




2 


k 


2 


+ — 1， 


: k 






2 


2 


r 々 


2 


J (4) 


m 




i ， 2 ^ = 0) 


l 


m 


「軏 


+i 


k f^t+i 


2 

々"■ Hl 

" T " 


+ 


~2 

m h 


— U 


Hh 




用 S 2 的方法，可以算出误差〜以及 
由 （1) 可以算出所有的 
令 


进而 


F 0 = F 


F 


m4*l 




■ fflt+l 


m 


F m + 内 i +1 F 

^ 0, 


C 5) 


-m -m ¥- 


则误差为 




F 


* 




读者可能会问，是否还有其他方法使得误差（或精密度) 
更小一些？这个问题的回答是否定的，对于单峰函数的全体 
而言，精度不能再好了.这里不再详细讨论，请读者自己证 
明，例如，可以参考第一章 S 5的方法，先建立一批不等式, 
然后用归纳法证明. 


5 6. 试验批数不定的情形 


有时，虽然每批试验的次数一定，但是试骑的批数并不能 


- 5 t . 
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预先给定，要求试验可以连续做下去，怎么办? 
先讨论每批做偶数 (2/) 个试验的情形. 


由 S 4的例，每批作2/个试验，共做4批，第一个试验点 
在 

A 從 2 (/ + I)— 1 — 1 

1 2(/ + 1)^ — 1 ' 

第二个栻验点在 


+ 


^4 +i 


2 (f + 1) 卜 1 

20 + 1> — l fc 


当々 一 C 0 时，二试点位置都趋 于土. 由此可见，这个方法 


当々 — 0 Q 时是完全失败的，也就是说，预先给定的批数4 
越大 5 两点越近，分辨起来越困难，正因为如此，我们建议用下 
面的方法来安排. 

先把试验范围等分为 （2/ + 1) 段，在个分点上做第 
一 批泣验，比较 结果， 留下较好的点*及其左右一段，然后把 
这二段都等分为0 + 1) 段，在分点处做第二批试验等等，下 
图表明了/ = 3的 情況： 



图20 

■ 

值得指出，根据实际情况，我们可以用混合的方法来安 
排，也就是先用 S 1的方法，做到某一预定的程度后，还可细 
分时，再用上述方法. 

再讨论每批做奇数（2< —》）个试验的情形. 

* 5$ » 
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由 S 3 的例 h 每枇作（2< — 1) 个试验，共做々批，第 
个试验点在 


这 1 


F 


k 


F 


4+1 


2 


+ 3 



4 


+ ^ t(t + 4) 


k 


2 


3 



4 



0 + 4)1 卜 1 、 


2 


2 



Z + V ’ i ( f + 4)1 

+ 4 / 

2 J 


—fl — 3 
2 \ 

当 o > 时,它趋于 



t(t -b 4 ) 

~2 



Q 


2 



2 



4 


t 




(/+ 4) 

同样可以得出，第二个试验点与第一个试验点 的距离 

F 


ai 


k+t 


当 CC 时 3 它趋于 

2 


衫 2 [ 


2 


+ V ’ 十 +) J 


2 


+ 2 



宁) • 


容妫看出 


6-{-0 2 


t 




这样，对每批做（以一1)个试验，我们可以按上式取 t 即 


6 


2 


— 1 + 



+ 4 


先把试验蒞围分成2/段，其长度比是 

行：沪 :(9: 屮 ： .. 

在分点做第一批试验后，余下是 
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e：e 2 或伊 : 6 . 


由亍对称关系，只讨论前者. 
因为 


t{e 2 + e 3 ) — te(e -h e 2 ) = e 9 

故把长为 0 的一段分为 2/ 段，其长度比是 

汗：少：铲：.…：史：护. 


再在分点处做第二批试验，余下是 

屮:史或； 6^:6\ 


因为 

K 仔 + 0 1 ) - /(a + 浐）汗 一 屮， 

又可 把艮为 屮的一段分为2/段，再做笫三批试验，等等.如 
此做下去，直到合乎要求为止. 

这个方法是黄金分割法的推广，/«1时，就是黄金割法 
(一次做一个试验)， 

最后，看看这个方法的精度如何？ 

记为每批做〖个试验，共做 K 批后余下的区间长度. 
对于 < =2/的情形，有 


n 2 ° 


2 


It 




2 


It 





2 


It 



对于彡 = 2/ 


1 的情形，有 
= & + d 2 


以卜” = 屮 + 屮 




0 
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= 屮 + 屮 





=妒 + e k+ ' -= 



§7. 是否最好的安排 


上面几节讲的分批试验法是黄金分割法和分数法的自然 
推广,因而第一章里的许多槪念和方法都可以平行移过来用. 

现在我们对 S 6的方法作进一步的讨论. 

仍然先讨论每批做^个试验的情况. 

我们的方法着 眼于： 做一批试验以后，进行比较,只留下 
一个区间，其中包含已经做过的、结果最好的那个试验点，也 
即如下图的形式.我们将它放大适当倍数，使留下的区间为 
(0, 1 ), 其中0 处是做过试验的点，现在的问 题是： 再同时 
做2/个试验，这2/个点布置在何处，使下一次留下的区间尽 
可能短？ 


图 21 

连同上次留下的已试验点，共有+ 1) 个试验点 ，设 
安排在 &2,…， 处，如果最好的在心，则下次 试验将 
在 

d ^ cece ^ 

中做，因此下次试验的区间长度应当是 

T 由 max (^h-i — % = 0 ， On+i ^ 1 . 
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由于 


T — max C^p+i — ^*-[) ^ rmx 0^ — 

l<i ? 或 Ji 1 十 1 l<r«^+l 疋 


> 


1 

— T 2 j 

"l 丄 Jr = 1 


〜 一 2) = 


(^+2 — B q ) 


也就是最短的因子是 


，比这个数更小的就无法保证，显 


然，均分法即 




1 


p 20 -h ir 

就可达到这一目的，逸就是 S 2 所说的方法， 

当然，如果第一批试验不是按等分取点，而是更细致些， 
安徘为 、 


f a 



b a f 

0 



L 


图22 

以后再桉等分取点.容易征明这是最好的安排，其中 

1 h 十1 


b 


2(7 + I) 1 - r 


a 


Y(T^ry~T 7 


满足 


( 礴 + b 、 i 七 a 


2 / + 2 


2( / + I ) 2 — 1 


•# 


1 


1 


1, 


2(/ + I) 2 — 1 

以后余下 

alb b'a 

把长为《的一段等分为 （2/ +〗）分，即每段投为 h 在分点 
上做试验，以后均按 S 2的方法.这样第一批试验后余下区 
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间长为 


l x = a ^ if 


2 / + 2 


2 (t + iy 一 1 


< 


2 


It 


以后缩短的因子都是击. 


为了取点简便，而精度又不相差太大，我们建议采用 s 6 
介绍的方法. 

再讨论每批做 （ 2 / — 1 ) 个试验的情形 3 
这时，可以证明 S 6 给出的方法是最好的.为此，我们只 
需证明下面两个引理，然后按照第一章 S 5 的思考方法去论 

证 ■ 

记 A?— 13 为々批试验后的精度，我们有 
引理 1 当免> 1 时，恒有 

证左迠不等式是显然的，故只证右方的不等式 . 

记 Od， 是第々批试验余下区间之最大者，即 

太糾 一〜 = 处卜 0 + ⑴ 

第（々+ 1 ) 批试验 （ 2 / — 1 ) 个，连上次留下的已试验点，共 
有 2 /个试验点，记为&，^ 21 ， 

且 

尤 0 < 尤 t C h C • • _ < x 2 t < h+1, 

我们有 

x^ x — x p ^ < (p i ， 2 J … ， 2/), 

由⑺ 

f - 1 

At ^ A 1?+1 一 A = ( X 2^hl — 尤 2r) 十 (Am — ^2t) 

f ■= 0 

< / W + (〜+! — 〜)- ⑴ 

在区间 0 2 „ r i /+ l ) 中，第 U + 2 ) 批试验点最多只能安排 
( 21 - 1 ) 个， 
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^+1 — ^ 沈 0 . ⑷ 

由 （3) 和 （4) 便得 （1) 式. 

引理2 当々会1时，有 

-h m (5) 

其中{巧巧由递归公式 

朽0=1， = - t ( F % + F ^ U ) 〔6) 

给出. 

由引理1，用归纳法易证. 

对于 S 6给出的安排方法，记 ❾ 为第4批试验后留下 
的区间投度，容易验 i 正 

没卜 0 - Klf ^ 4-财 ， Oh ⑺ 

" =，(平從 + a > Un > l ). (8) 

根据（1)，00,（7)， （8), 用第一章 S 5的方法易证， 
定理对任一每批做 一 1) 个试验的安排方法 
记 Ar i 为第々批试验后的精度，若^不同于 S 2的方法 5 则 
一定有，，使当尺 >时， h 

△浐™ 1 ) > ， f _1 ) 口 

即： 5 2的方法是最优的安排方法. 

5 8. 依某种要求进行试验 


每批同时做几个试验虽然有加快速度的好处，但是也有 
增加试验费用的坏处，并且每批试验个数增多，原材料消耗 
比按正比增加得还要快.例如，用分数法89/144 做 十次试 
验的效果和一批同时做 H 3 个试验的效果相等 5 时间缩短为 
1/10,但试验的代价却是 14 J 倍，因此，必须权衡利弊来决定 

j 
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采取何种方案.此外，一批做多个试验时，试验范围缩小得很 
快，例如，每批做7个试验 3 共做六批，精度为2/19776,即为 
万分之一.如上所述，两个试验点愈靠近，则通过试验来分辨 
它们也就愈困难，因此考虑批数做得很多，意义不大.遇到这 
种情况，首先看看试验区间要求分得多么细，再决定一批做几 
个试验，共做几批，来安排试验点等等 . 

例如，有八套仪器 7 可以同时做试验，但在试验点比较接 
近时，要泶重复试验来避免误差，简单地说，我们希望这样安 
排： 进行四批试验第一批八个试验点在八个不同的地方做， 
第二批四个试验点在四个不同的地方各做两个试验，第三批 
两个试验点在两个不同的地方各做四个试验，第四枇在一点 
做六个试验，并在留下的较好点补做二个试验加以比较.也就 
是说每枇试验次数依次为（8, 4, 2: 1). 

如果按 S 5的理论上推导，容易写出 

广 1 ^ 

at = 3^2 + 2^4 

= 2 a ^ -h ^ 

I ^3 = ^4 + 

议 4 tf== 

并算出误差 A = l / S 9. 但是，这样就要求第一批试验点两两 
相距1/89,与我们的前提“试验点比较接近吋，要求重复试验 
来避免误差”不相符合了，我们可以用下面的安排 方法： 

0 先分试验范围为九等分，在5点上做八个试验（试验 
点相距 1/9) 

人 1 it. > > _ I r I I 

U i 


2) 对余下的区间仍 按等分 (试验 点相距 1/27) 
r * 


PDF 文件使用 11 p df F ac1 0 r y 11 试用版 本创建 www. f i nepr i nt. com, cn 























每个打 * 处各做二个试验. 

3) 对余下的区间分成六等分，在第2, 5分点各做四个 
试验： 


4) 最后在余下的一点做试验. 

这样安排,误差是1/9 X 1/3 X 1/3 = 1/81,与前面安 
排的误差 1/S9 相差并不大，但却能有效地避免了两点相距很 
近不易分辨试验结果的困难. 

§9. 重复性试验的分辨问题 


我们在上一节说过，有时需要用重复试验来避免误差，由 
于在一个点通过试验所得出的 / GO 的数据，并不— 
定是 /OX 而是可能有误差的、例如测得的 y=/(a) + S， 
?就是误差，为避免误差的严重影响，我们往往在同一条件 
X = a 的情况下，多做些试验，得出以下的数据 

yu yz * * mw ^ 

在另一点 x = b , 也做了些试验，各得出 

如果^都比 〆大，则两个试验是可以分辨的，也就是可以指 
望/0) > fQO - 但是，有些 h 比有些/大，而有些 h 比有些 
M 小，怎么办？ 

第一个建议是用平均数法，如前，在^ 时所测出的数 

据的平均数 
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歹 = 丄 （)、+ y 2 + + …+ > V ) 

y f ^ — ( y ; + y ; + ■… + y :) 

大，就认为 /(«) 比/幼）大， 

第二个建议，看是否是 

y > 每一个 / (卜 1，2，… 0) 
7 < 每 一 个 (/ 1 i 2, ■ * w ), 

如果是这样，我们认为 /( a ) > Kb) 的把握更大了. 

第三个建 议： 考虑不适合不等式 （1) 的个数如果 
中有 f 个不适合（0,则我们可以闱 

1 P 


来检定 1(a) > fOO. 

第四个 建议： （ h ， 〆 ） 共有^对 ，其中 适合于 

yi > y ) 

的有0对，则/00 > iih ) 的把握是 

St * 

例如在两点各做了四个试验，得出数据 


yi ^ 

- 1.19^ 

^2 == 

i.2i ， >3 = 

= 1,21 ， 

少 4 = 

n 


- 1.10, 

)4 _ 

1.15 a >3 M 

=1.18, 

yi = 

1.20, 


我们只要依大小排队 

y\ < < yi < ^ 

h 与 〆 .共有 4 x 4 = 16 种 关系，15个都适合 〆 < h ， 只有 
一个 a < y ；, 因此 Kdb ) 的可靠性是15/16二 0.9375 + 
上述这些处理方法，是比较简便而又避免了不必要的繁 
琰 的处理方法，是比校容易掌握的.但在同一点反复 做试验 
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得到数据很多的吋候，我们不排斤运用“正态分 布” （如果数据 
大致服从此规律时）中的标准离差等方法来判断最优点出现 
的可能性的范围. 

还有一点值得注意，接近于 m . 难以分辨吋，即使 
就是时（或差不多相等）， /0 O 的最大值在 o , 
O 之间的可能性也是很大的，换言之，〃之左， i 之右都可以 
舍弃，因而建议，暂且不管之外，只对 M 开始用 
黄金分割法或其他方法 5 这样可以省做不少试验/ 


S 10+ _ 卩单峰 的情况如何办 

I 

■ ■ 

有人 会问， 前面所讲的只适合于“单峰”的情况，多峰(即 
有几个点，其附近都比它差)的情况怎么办？我们 建议： 

K 先不管它是单峰还是多峄，就按单峰的方法去做.'找 
到一个“峰3后，如果符合要求，就先照它生产，然后有时间继 
续再寻找其它可能的更高的 “峰” (即分区 寻找乂 

2 . 先做一批分布得比较均匀疏属的试验，看其是否 有“多 
峰”的现象出现，如 果有， 则分区寻找.这时 5 第一批试验的点 

最好依以下的比例划分 

<x：^ = 0 . 618:0382 

ce _卢 a a fl 1 A a _ _ 


图 23 

例如，三个分点，可以取之如 

I_ a 〜 乒 * a ^ 


图汾 

^ 69 m 
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使留下的成为 

-^- 0 _ i _I 或 * 一《 £ — _ o- —— - * 

m 25 

的形式，接下去便可用 0.618 法了. 

最后要说明一下，有些问题，并不是寻找单峰函数的最大 
值，而是如对分法那样，寻找一个合适的点.在区间 04) 中 
有一个 单谪增 加函数 /0 O ，还有一个数 G 

/(5) < ^ < / O ). 

要寻找 q 使 

/(4) w 入 

我们知道，对分法是在区间的中点做试验，若 f 

则舍弃若则舍弃 

再在余下区间的中点做拭验，如此继续做几次试验后，余下区 
间的长度是 

l n 

容易证明，对于全体单峰函数来说，这个办法是最好的.如果 
条件 允谇每 批做几个试验，那又如何安排呢？很简单，不管 m 
是奇数还是偶数，这时都是用均分法最好，例如每批做三个实 
验,就把 ( ad ) 四等分， 

I - ^ -- >f- — at — 1 

a &t di 轟 

图26 

在三个分点 A ， 仏 处做,然后再考察试验的结果，例如 

t 70 • 


PDF 文件使用 11 p df F ac1 0 r y 11 试用版 本创建 www. f i nepr i nt. com, cn 

















f ⑻: > c ， 则舍弃 （《 ，民） 和（馬，办），再四等分 
(‘色），在分点做第二批试验，如此继续做》批试验，余下的 
含有 A 的区间长度是 * 

士 0 — a). 

4 n 

也容易证明这个安排方法是最好的. 
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第二部分多因素优选法 


这一部分考虑的是多因素优选法，我彳门将分章研究这些 
方法在多因素情况下的-广.和前面一样，每章都有一段方 
法的述要，在明确了方法的使用以后，再逐歩深入研究方法的 
收敛性和一些有关的数学问题，最后一章是讨论优选法与计 
箅数学的关系.所介绍的优选法，实质上就是求最大值(或最 
小值）的快速计算法.不仅如此，很多计算问题都可以化为极 
值问题来解决，而极值问题离散化后，就与优选法有关. 因而 
优选法在计算数学上也是很重要的. 
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第一章双因素优选法 


k 



述 要 

绝大部分的双因素优选法都可以推广到两个以上因素的 
情况，只是所需的工作量> 将随着因素个数的增加而迅速增 
长.下面我们先讲 方法： 

1) 网格法把因素 t 和 y 的优选范围（区域） 

a x <Z c <Z y <1 d 

各分为【和％个等分，于是总共有^个格子，在每个格子中 
做一个试验，共做^个试验后，哪个好用哪个，这个办法很 
简单，但它的工作量比起一个因素来说，是大大增加了.例如 



10,就要做 1 C 0 个试验 


图 27 

假定处理的是单峰问题，也就是把^ y 平面作力水平 
面，试验出来的数值 s 看成为这一点的高度.这样的图形就 
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图 30 

像一座高山，我们假定这一片地区只有一个制高点，在平面上 
画出等高线后，这些等高线形成一圈圈的互不相交的封闭曲 
线.我们还假定，每一圈都是凸的，就和树木的年轮一样. 

2) 对开法裁一块方格纸代表优选范围 

a <x < if/ c Cy <d w 

左右对折起来，在中线 

x Hh 方) 

2 

上，用单因素法找最大值，设在 r 点取这最大值，再上下对 
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折，在中线 


y — 丄 (V +乃 

Z 

上，用单因素法找最大值，假定在&点取这值. 

比较 P 点和0点的结果，如果 P 大，裁掉纸的左半张(不 
然裁掉下半张 )， 再用同样的方法来处理余下的半张纸，不断 
地去其一半，逐步迆得出所需要的结果4 

3) 旋升法先在一条中线（例如左态对折线）上求最大 
值，假定在 A 点取这最大值，逋过 A 作水平线，在这水平线 

- 75 • 


PDF 文件使用 11 p df F ac 1 0 r y 11 试用版 本创建 www . f i nepr i nt . com, cn 








































上找最火值.假定在 A 处取，又在逋过朽的垂线上找最大 
值，而在巧取这最大值，等等.继续做下去，存荫影的部分以 
后不 要再考虑了. 

4) 平行线法两个因素中， 一 个例如 x 易于调整，另一 
个例 如 > 不易调整，则建议用“平行线法 '先把 y 固定在范围 
心， V )的 0.61 S 处，即取 

y = c — r)0.618 ， 

用单因素法找最大佶，假定在 P 点取这值，再把 y 固定在范围 
0, o 的 0*382 处，即取 

y — c (d — r)0.382 ， 

用笮因素法找最大值，假 
定在 P 点取这值，比较心 
2的结果，如果 P 大，则去 

掉下面部分，即去掉 
y < c + Cd — ^)0,382 

的部分（否则去掉 h 面的 
部分），再用同样方法处理 
余下的部分，如此继续.这 
个方法也适用于分批试验 
的情况. 

5) 按格上升法把所考虑的区域打上格子，与2)，3乂 
4) 类似，但用分数法来代替黄金分割法. 

6) 翻筋斗法从一个等边三角形出发(如图 36.), 
在顶点各做一个试验，如果在 C 点所做的试验最好，则作 C 点 
的对顶同样大的等边三角形 CDE . 在 Z ), £处做试验，如果 
D 点最好，则再作 D 点的对顶同样大的等边三角形，…，一直 
做下去.如果在 F ， G 处做试验，都没有点好，则取 FZ ) 及 
GD 的中点 f 、以做试验，也可以取 CD 及的中点做试 
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验，再甩以上的方法，如果在0的两边一分再分都没有 D 点 
好，一般说来， D 就是最好点了. 



A 



§1•对开法 

实际上，我们的优选问题就是在区域 

a x <i b y c <i y <d 

内求 ar == /O, >0 的最大值，所谓等高线，就是固定〜由方 
程 z /O, y) 在 Or, 妁平面上所定义的曲线.我们假定等 
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髙线有连续切线. 

在 )0 平面上，一条直线上取最大値4的点，一定是 
等髙线 Kx ， y ) — e 与这直线的切点，因为这直线不可能与 
它再交于另一点.由于等髙线是凸的，所以，使 a 的点一 
定分布在这直线的一侧，现在用这个原理来观察对开法. 

(1) 这方法一定可以逐步通过最大值 C 或简称收敛法）. 

我们看到，应用这个方法 

做了第二枇试验后，每做一批， 
区域的面积减少一半，因此收 
敛性是没有问题的，但有一个 
特殊情况要注意，就是会不会 
出现 P 点和0点恰巧都是中点 
的 情况， 即橫线纵线上的值都 
没有中心点的值大，这说明如 
果等高线 

Kh )0 隠 d 

不缩成一点(制髙点），则等高线在这点不存在切线，或者说切 
线在这点不连续. 

(2) 有效性 

如果单因素是用黄金分割法进行的，假定所要求的 
精密度各为 e 及^经造当的放大缩小，不妨假定在 
第一条线上作/次试验后，所达到的精密度是 

一 心如卜 1 ， w — — — 1 + 5 ) , 

- - 2 


p 

1 



[ Q 


图37 


现在所要求的是 

C 々一 a ) 如卜 2 > S 会（方 一 U fl ~ l J 

即 
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log 


b — 


£ 


> / — 2 , 


log 


也就是；差不多等于 

b 一 


log 


s 


V ]os i 


log 


+/ 1 嗯去 • 


+ 


以后每做一批，区域®积去掉1/2,经 f 次单变量处理 
后，得出的面积是原来的 


* 


2 P 


如果做到 


一 一 < e 2 < g — oty — o 

2^ 2户— 2 

当然就满足要求了，也就是 p 大致等于 


2 log 


8 


tog2 


因此， 


2 


log /(】叫 


log 2) 


次试验可迖目的，也就是当 e 很小时，无穷大的阶是 


log 


2 


£ 




比起网格法要好得多. 


e 


附记在 P ， P 二点做出来的试验数锯如果相等（或无法 
辨认好坏），这是件好事，因为这说明 P 点和 P 点位于同一条 
等高线上了（见图38 ) ，所以可以把图上的下半块和左半块都 
去掉，仅留下第一象限.这一点对统计方法来说也是十分有利 


* 
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的，当两点试验数据的可分辨性十分接近时，我们干脆丢掉面 
积的 3/4. 



在进行旋升法的时候，请 注意： 划掉一个区域的原则，也 
就是有了 A 点与 Pz 点之后，划掉通过 A 点的直线所分开的 
不包含巧点的部分，这样可以少做不少试验，当然， z 值不断 
上升的情况可 IS 是外向的，也可能是内向的，也可能兼而有之 
(图 40). 

从方法来看，每次求最大值，可能比財开法好些，但并不 
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能保证每次划去的邵分大宁或等宁原籴面积的一半.也就是 
说，上算的可能性虽然很大，但不上算的可能性也是不能排除 
的，不象对开法稳拿一半. 

(1) 收敛性 
如果一直做下去，则 

Kpi ) < & = KpO = /(?„) 〈… 

是一个单调上升的数列 ，而 由问题性质看出这数列是有界的， 
因此 

lim ^ 

我们画下高度等于 

…… 

的等高线 •* 、如果 S = %不是制高点的高 
程，则有一条等高线其高程是当》充分大吋，^在 
A 的法线将穿人^之内，换言之，比 A 更大，这是矛盾. 
朽在 Ow 之中，所以趋于制高点. 

C2) 有效性 

要使误荖小于00) ? 需要的试验次数的数量级是 



证明请看下一章的 S 1. 

S 3 .平行线法 

从对开法的讨论中立刻明白平行线法一定逐步逼近最大 
暄.如果单 H 素是用黄金分割法，则每条直线上要试验的次 
数大约是 
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而每做一条线上的试验启 ，区 域面积缩小沟原束的这样> 
次后，留下的面积是原来的 

当要求 * 和 y 都达到 s 的精度，的阶数约为 


因此一共要做的试验次数约为 

( iog t)V ( iog iT- 

同样可以讨论分批试验的情况. 

■ 

S 4, 两个因素的离散情形 



例如优选的范围是一个21 X 13的格子，我们现在不用 
紂开法，因为平分下架不一定在格子点上，也不一定上算，怎 
么办呢？我们在 r =13的直线上用分数法做五次试验.又 
在7 8的直线上也用分数法，这时 T 点已做过试验 ，因 此做 

6—1 次拭验.各得一最优点 S 记之为 P 点， 0点.比较 P 点与 
1?点，如果 e 点比 f 点好，则裁取一个3 X 13的格子（图42乂 
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在这个图形中 ，取^ « 13 + 5 = 18 还是取 r — 21 —5 - 16 3 
可以这样来 考虑: 因为 1 S 更靠近13点些，因而取 
18, 再在这条线上用分数法.这样做并不能排除不上算的可 
能，但上算的可能性大些(因为近“大”处得“更 大，. 

如果在每个格子点上作试验，共 
需试验 20X12=- 240 次，闬现在的 

方法，最多只要30次就可以了，如果 
纵横格子个数并不恰巧等于某一 
那末可以添上一些或冒点风险减少一 

些，以凑成例如在 

0 <x < 1 S 

吋，不妨添上些成 D < x <2 i ， 或减 
掉些，而成 

0 < jc < 13 , 



图 43 


当然也可以用平行线法，读者自行补充 . 


S 5.翻筋斗法 


当因素的变化不能跳跃式调整时，此方法有其优点 ，翻筋 
斗法的理论根 据是： 通过一点的两个方向的微分都等于零， 
一 般来说，这点就是极值点. 

其实前面关于等边三角形的哏制不是必需的，根据具体 
情况用直角三角肜，任意三角形都行.如果用四边形， 五边 
形，…来翻筋斗效果如何呢？例如，有人建议用矩形，在矩形 
的四个顶点做试验，然后依最好点翻过去，在新的三个顶点 
做试验，等等，称为“矩形调优法' 这 样做好不好？这样既多 
做了试验，又没有很好的理论根据.因为这里同时考虑了三 
个方向，而极值问题一般只要考虑两个方向就够了，因此“矩 
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形调优 h 就不如用三角形翻筋斗 i 翻筋斗法不难推广到多个 
因素的倩形，例 如:三 个因素，用正匹面体在四个顶点做试验， 
如图所示，依最好的一点翻过去， 等等. 



图妇 

关于” 个因萦则是找出一个〃维的正单纯形，在它的 
«+1个顶点做试验，怎样找顶点？我们建议用以下的 方法: 

设3点的坐标是0^, 以 j 为顶点边长为 A 

的”个顶点 A ly A 2 , 的坐标可用下法得之，先做月 

点： 

b t = -f- ah ^ bi = a 2 ^r ah^ * ■ * 5 if n ^ -h ah^ 

这儿 




也就是每一坐标上都加上 d ， A 就是 S 的第一个坐标加上 
4二，冷就是5的第二个坐标加上一儿就是 S 的 

V 2 V Z 

第《个坐标加上其余坐标不变.证明是不难的^ 3与 

V 2 


4的距离平方 等于： 

(?j — 1)( 。 A) 2 + fdA -h — ^一 



一 i) 找 2 + W + ^/ 2 a -t — 

2」 
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\/ 2 a + 



(\/ n 


1> 


n 




〒 h 2 , 

所以 d 与汰的距离等于 A. 
又冷与 AiH ^ 0 的距离等于 



另一 方法： 假设^是原点， A « 1 (都不失其普遍性），取 

山=( 1 ， 0 ， 0 ， ’-. jO )， 

■ 

A 2 = ( Aj … J 0)， 

由于 a 与^的距离等于1， a 与爲的距离也等: n ， 所以 

朽 + 4 = 1 3 (ft — I ) 2 + ^2^1. 



a 2 ^ 士^取 ]£ 号得 

2 



我们可以 一 维一维地添加 ，一 般有 

Ai *= (ft，A，* * * s 參 «；s 0, . * ■，0)， 

这儿 

- A ~ iM V^crro 5 ai = ^i l ir' 

证明： d 与為的距离平方等于 

珩 + 的 + + * . 十 fiU -h a ) 


1 1 十 

i -h i 

i 卸 （ ？ + i) 

21 


S5 
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1( 告 


2(/? 十 O 





又4与 ASI > 的距离平方等于 

C 〜一卢™) 1 + 券 h + •. * + / S| Ml + a] 







2 pCp 1) 




m 

2 (m 4 - 1) 



2( 芦 + 1) 





2(m + 1) 




2 (^ + 1 ) 



+ — + — = 1 
21 2 

即得所证， - ^ 

新的试验点是很容易求 出的： 在 （《 + 1) 个顶点杓， 
p 3 , •…， 上做试验，其中 

Pi =" C ^ f lS ^ J l5 …，〜） G b 1 J 2，…， n + 1). 

假设巧上的试验结杲最好 ，依巧 找出各顶点的对称点 

«1,2 ? - ■*,« + 1)： 

即 

^ — 2P S — Pi , ( j . = 1，2 ? …+ 1乂 

然后在新的顶点钺试验，找出最好的 一点， 又依这点找对称 
点 3 —直做下去，如果做不下去了，比如 A 仍是最好点，这时 
将距离缩小 3 例如缩小 一半： 

H 
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^ - p } 




(d) 


即 




在 p ；% &，•••， 继续做试验 i 


■ 


窜 7 


- 
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第二章最陡上升法 


述 要 




先在 + O 个点 〜 •. ■〜）， （☆ + 5 1S 

尤"）， * * -(^ 1? x l 7 …，丈厂 1? ar „ + <5 n ) 上各做一次试验，得出 


数据 


苽 I, 怎 


⑴ 


Z 


在线段 


⑴ 


X 








s n 

上用单变数方法找出 x 的最大值，这儿 * 的变化范围由以下 
的不等式决定 

C0 ^>0, 


U ) 


(ii) Si ^ Xf -h — O 




,2 




* w * 




如果找到的最大值在 f & 时取，则令 






以 （ h ， "•，代替 …， 继续进行 ，方的值 
不断增大，一宣做到 （ h ， r 2 ，• * _，与 （: Vo A ，' ' *> )0 的 
距离小于预给的精密度为止. 


§1. 最陡上升法 

假定在区域« < x ^ ^ j 内， 

* s $ - 
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z 


= F <： x , y -) ⑴ 

是一个向上凸的曲面，处处有二阶连续微商，并且假定在讨论 
的区域内仅有唯一的制高点. 


在制高点 Ofl ， 外）处 3 


^ ^ y ) 


C2) 


一 定适'合 


尸 * Oo , yo ) = 0， F y ( x a , y a ) = 0 # 

由于曲面向上凸，所以 


(3) 


&2F 阳：^ 


a 2 F 




2 


dx 2 dxdy … dy 2 

是一个定负的二次型，也就是一定存在两个正数 A / 与 m 使 


Ox 1 






dy 


2 


dxdy 

< 一 wx 

由闭区域内连续函数的性质， A / ， fn 与 X ， 无关1 
从一点 Ca ，）0 出发，若 F 0 1Jf >o = a 对等高线 


⑷ 


FO, y) 


方1 


(5) 


沿着法线方向作直线 


J ： * JCj 




y ^ 


ay 


0 ^ OX (6) 


用单变数法，沿此直线方向找 


£■ W ^ t 


idf ) 

Va 7 


yi + 


0 F 

dy 



的最大值，假定在 


h 吋取最大值，令 


方 2 


gOd =尸 (A 


\ 


dE 

dx 


w+dfW ⑺ 


F ( x 2 y yi}y 


由于 K *) 也是凸函数，因此 


* 69 


晷 
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d 


dt 


^0) 


d _ 

di 


F U 


OF 

a 7 


yi 




o 


■ 


CO 


我们现在硏究(与制高点高程差的进展比例) 


z 


0 




土 — 尹 （¥ ;Vq) ― F0 2 , y 2 ) 


H — 


FO 。， y c ) — fO"jO 

尸0^， y ?) — FCu ) 

■ ■ 了 ■ 

尸 ( A ， 外 ）— 


(9) 


1) 分母 


F ( h ， y 0 ) — F 0 1? v L ) 




Jo de 


+ i ( x { — r 0 ) 


y 


Kyi — 





dt 


d 2 




0 ds 


2 


FCx'a + — Jr 0 ) 


+ jr(yi — yo)ih 




- i : 


dt \ [C^j — jf a ) 2 F 


X 


2(h — X 0 )(y y — y a )F” 屮 — y^yF/^s^ 

由 （4.) 可知 

FG. 。， yo) — ^'C^i 3 yi) ^ MIX 尤 i 一 x 0 ) 2 




Oi — ya) 2 ] 1 dt 

0 


ds 


o 


M 


* — IC^i — ^a) 2 + Cvi — yo) z l. 


( 10 ) 


2) 再有(根据 (7)) 


FCh ， y 2 ) — FCu) 


» dt 


F x 


f )/ 




yi 




dF ^ 

Qy 

dF _ 

dy 


\ 






dt 


d 2 


0 


ds 2 


F x 


dF 

dx 


i 



^ =~= — 


dt 


r (( 


dF 

dx 


2 


F 


x 


_ 


90 
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2 


(-t)m 


F 


m 






m 


= 


m 


2 


( f ):] i : 

[(H 


Oi — r)dt 


3 ) 从 （ s ) 式推得 


if)：]- 


00 






h + 


pF 

d7^ tly h 


我们把 


(fMm 

( if )/ 


F 


y 





dF 

Qy 


02 ) 


F 


X V x i 


OF 
dx 


A ， A + 


OF 


dv 


d 

dt 


F 




i 


及 


dF 

dx 


jp 


yi + 


匕(: JJ5 ： Ki ) 

(f A , 




® J - 


o df 


F 


y 


A + 


dF 

dx 




y\ 


代人 C 〗2 ) 可知 

咖 m-m 


dF 

dy 




dx 


d 


F 


( f ) 


a 无 


^yi + 


pF 

dx / 


> p i + 


P.^\fdF 


dx 


dy 


i 


^)M (f), 
德 ) i . 

m 


d 

•" !■ 

dt 


F 


y 


x 


w 


0 



dF 

dx 


F 


X 


£ 


F 




F 




1 心， 


03) 
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由（ 4 )可知 


m 


dF 


2 


dx A 


<A/ 


f 

、. 

OF 

Qx 


dF 

Oy 


21 


1」 


^ (S ^ 


dF\ 2 

dy 


2 


dF 

W 


i 


尤 I ， 


即得 




_ 


( 14 ) 


4) 由 （ 3 ) 可知 


dF 

aT 


Kh yi)— 方 o, yn) 


1 d 


o dr 


— A)， 九 + ^C} ; i _ y^))dt 


及 

dF 
dy 

并推出 


1 


1 d_ 
0 di 


Py(^a + 心 l — A)，Yo + Ky — )0, 


dF 

dx 


(^i — ^0 


dF 


Cyi —妁） 


由于 


JT 


3 


■- dy / 1 

+ 2( 太】一太 0〕C，1 — >0尸 jry 

Cyi — y o yF^}dt 
2 ^ F, y 4 -铲~是定负的，所以 




dF 

dx J 


O — Xo) 


dF 


C^i — ya) 


0 


I a 


〜) 2 iV 十 2C^i — hXyi — 7o)hy 
y^yFy^l dt > — x o y + Cyi 




Oi 

y0 2 L (15) 


利用 Schwas 不等式得到 


一 A) 2 + Cyi — jo) 2 ] ^ 


dF \ 2 

d 7 /i 
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1 


dF \ 2 ^ 

dy n. 


[O 一 + Oi 一 


于是 


抱 2 [( 方 1 — a) 2 + Ca — yo) 2 l ^ 


dF 


2 


dF 


2 


06 ) 


综合 1)，2), 3) 及 4) 的结果，可知 


FCjtz, y 2 ) — FCx^Vi) 


> 


2 


mt 


2 


8 F\i 


dx 


dF 

dy 


尸（龙 a ， — FOu yO i 


2 


m Kx^ — x 9 y + c^i — ^o) 2 I 


/ 0F\2 


> 


m 


\ 


dx 


dF 

W 


m 3 c^i — x Q y + c^i — y&y 


> 


m 


3 


M 


也就是 




您 i 


< 


m \ 


Zft _ z 


i 


M 


/• 


(⑺ 


这个结果值得注意的地 方是： 右边与％无关，也就是从任 


点开始，做一次后，与制高点的高程差缩小1 一 
果做了 /次，则高程差与原来的高程差的氏例是 


m 


M 


i 


倍，如 


(-©7* 


取/使 


©1 




tn 


M 


3 


则必 


_ 
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也就是 / 的数量斷 


在一条直线上做黄金分割法的数量级也是 log+， 


而每 


找一陡升方向需要做二个试验，因此总共的试验次数的数量 


级是 



这方法可以推广到 f 个因素，而由于每找一次陡升方向 
需做/ 个试验，所以总的试验次数的数量级仍是 

( Io 8 t)- 

最陡上升法形式上与第一章所讲的旋升法有昕不同，实质上 
是一祥的.旋升法的要点是先求一条直线厂上 4 W 的最 
大值，敢这最大値的点就是这直线与某等高线的切 
点，陡升方向就是法线方向，也就是在 Otu yO 点找过这点的 
等高线的法线方向。因而原则上并无不同. 

旋升法不但易行，而且不必要用差分逼近微分，不必要每 
次做二个试验来确定陡升方向. 


§2.渐近陡升法 

在上节所讨论的最陡上升法中，如果函数 Fix . y ') 的表 
达武已知，二个偏导数可以通过计葬直接得到.但 

ox oy 

在实际应用时，对所求的目标函数絶大部分是连续函数的形 
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式， 是不知道的.因此，要求出在一点的陡升方向也只有用做 
试验的办法，一般可以这样做：在点 O , y ) 附近，沿坐标轴的 
方向各找一点 O + ）0 和 O , y + &)，在这二点做试验 3 

并计箅出 

f (， + 在 ” y) — f(x> y + 衣 j) 一 f(x 3 y} 

A 匕 1 ' 

当 & 取得适当小时，它们便给出点 O , y ) 处的两个偏导 
数的 C 近似〕值. 

把这个方法过渡到一般的情况，便是衡近陡升法了. 


保持本章开姶时述要部分所用的记号. 

我们现在每次沿最陡的方向上升距禽为 s 的一步，也就 
是取 〃使 


尤1 





与 Op …的距离等于 


0 > 0 ) 



取 s 小到可以分辨的情况，则得渐近陡升法，这就是在山坡 
上，每次向最陡的方向走上长度为 e 的一步，这个办法虽然不 
如前法快——一步跨到山脊——伹是对于不能跳跃做试验的 
情况比较合适. 

一般讲来 7 我们可以取合适的单位使 


并且不妨假定也就是用 




A + 


Z 


O ) 






77一 一 ％) 2 + 


e 


■ ■ 4 


+ (十一 o * 


(；™ 1 ^ 2 5 …， n ) 


来代替原来的 （ v i ， •*•，〜) 


95 
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& 3 .二 次模拟 


在 S 1中我们找到收敛因子 



在统计学中，我们经常(不太令人放心地)假定数学模型是二次 
的-下面将证明，如果真是二次的，收敛因子还可以有所改进，这 
说明在靠近最优点时。也就是用 Taylor 级数的二次及二次以 
下的项来逼近时，实际上收敛速度要比 S 〗所讲的要快得多 4 
也就是我们检査最陡上升法用在以下特例上的收敛情 

U ： 

n n 

卩 (*1，…， O _ ( + 2 乏 ] biXi — ^anXjXii 

这儿叫=叫，（化）是定正的对称方阵，它^特征拫 

0 ■< Zj ^ Ai ^ ^ 

经过刚体运动，以及在 F 上减一常数，我们可以不失普遍 
性地考虑以下的特例： 

p 

ft 

QCx ls …， o = — S “】， 

— I 

它的最大值是0,在 A = • * • = = 0处取. 

从一定点 Ot 严， d 〜）出发，在这点的陡升方向是 

C — 严，*"，一 方? )）. 

研究 

Hi ) ™ - 2 以％ •■•，：?)一 2 i rt x ?)0 

的最大值，由于 

/0) — 2 ^ 2 总 o a = - S 又乂°’ 

v =1 v =^X 

n n n 

+ 4/ S - 4 i 2 2^ - S 

tr m-X w =■! ¥ *1 


t 96 ^ 
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(s ^ 


2 


Z 

w =1 



E 奶 1 丫 
s 办 


所以，当 


时，即当 




时，取最大值 





VaV 妒 

j£ J /w-yrf* 2| 

vlvu ? 

/ J n/yA^i 




Q(4°, 4 1} ) — toy = — 2 i i?x m^ 


( s 奶: 


2 


ft. 


m 


I ] 办 r 

V =±1 

现在利用不等式（其证明在 下节) 


念受自 ^ < j - Q 


I 


X 


n 



I 

I 


2 


d 


得 


(±^)\ 

=1 / 


PC^ B .-.，a 





2 




A 


S ㈣ 



f: 



* ?7 
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回到一般的形式， 4 就是 M ，1就是 m , 因此收敛因子是 


M 


m 


m 




M — 


M 


m 



< 


(S 



由于下节末段所指出的论点，这因子是不能再改进的.这 
结果指出 si 的收敛因子是有改进的可能性的，如果能找到 
一个坫果连偏微商的差分逼近的误差也包括在内，那就更好. 


4. 反向 Schwarz 不等式 


现在证明上节的不 等式： 如果 


则 




~ < 








4 \M X n ^ X 

证： 此不等式的左边是 a , 的函数 


訄4 


^ p -b qX 


r 




其二阶微商大于零,因此在中不可能有最 X 值, 
仅当 A 或 h 时，它才能最大，也就是. 


-ti¬ 


ff 


^i a i 7 1 ^ A + 又 in 2j 


a t 


f 畔 


费 i 


71 


A-t i =?frt+i 


(IH 


7f 


s ^ 


I 


l 




Y) -S a * a 

*1 / i -f*+l 
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当且仅当 

■ 


及 

又1 = . ■ = … 

= 又, 

吋取等号. 

m fl 

2 a i 

/ i I i =3 W+X 






§又 收敛闶子的进一步改进 
在 s 1中我们粗略地估计出了 


< i - 

A — 芯 i 



⑴ 


而在 S 2中，我们对二次目标函数证明了：收敛因子1 一 

o /町 可以改进为在本节中我们将对—般的 

目标函数证明与二次目标函数相同的 结果. 

考虑 


怎 =/ 0 * 0， （2 ) 

这里 * = X ,， ... ，&乂假设 / o ) 具有三阶连续偏微商， 

且是一个一致的凹函数(在一维和二维时，即以前讲的“问上 
凸” 函数乂 

引进 记号： 


- 99 • 
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= v/Co = (H_L， …， 鲁_) (3) 

\ OX t dx 2 dXj / 

称为 / 在点 T 处的梯度 + 2 

HO) = ▽办）一 (4〕 

、 dxi dx / Kt ^ 

称为 f 在点 x 的二阶微分矩阵，又称 Hessian 阵. 

最陡上升法是从一点« O 产，+*•，亡 〕 ）出发沿着梯 
度方向作直线 

x — x ir} + tg v 0 > 0) (5) 

Cl — £(/—))，用单变数法求出 /O) 在这直线上的最大值. 
假定在/ =匕时取此值，令 

A 十 !>_〆 ，） + tvg ^ C6) 

然后从出发重复上面的步骤迭代进行. 

用51同样的方法，可以估岀 

fC^O - g ( x (0 > 4- /O ⑸一 ⑷ 

一 x rny dt _ J [ 发 O ⑴ + tCx^ - r co 0) 

一 g (，) K〆 一 x^ydt « ^ dt < y ，） 

- X ㈨ ) f O ⑼一 / 。 >))0㈤ — ，)) '心 

~ a i : ^ ~ 〆”〕("*+oc I -， I )) o ⑺ 

— ，)）' 厶 _ 丄（方⑻ — x ^) hx ^ v} — 

•2 

+ 0 (卜⑻一沪 P 乂 (7) 

这里札=//(户 〕 ），，是==/0)取最大值的点，故有 
K，）一0 ‘ 同样可得 

g , « O ⑷一 + fO ⑷一 x <9 ^dt 

= O ⑺一 x^)H v + 0(I〆 ，） 一 X ⑼ I 2 ) • (8) 
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所以 






因此 


含邮 + 0( i 〜呷) . 


又有 


Kx tv) ) — /0 (H ^) = — ⑻. 


gM(x iv} 4 - sg^glds 


=j v j ； 

■ 

=dt g v [ H y + oQg v \^ g^s 

=tlgMvgv + oCUvlO- 
2 

由于 Kx ^ + 在 f G 取极值，所以 


o =|/ C ^> + ^)I 


V 


g(x (p) + t v g v )gl 


+ o (\ g y \ 2 ')] S v y 


即 


t ， 成一 


1 > 




「 + 1心1) 


4 


代人 （10 乂得 


/ O ㈨ ）一 / O ㈣ ) 

綜合 （9) 和 （13)， 

/ C^ 0 O - Kx <^ _ 

/ O ⑺）一 Kx ^ 


1 U.l 


2 gjivgv 


/(/ ㈣ ）一 / O ㈦ ） 

- Kx^y 


IgJVg^^v + o (\ gA ^ 


g^7 l g v + oQx^-x^l') 


kJ 


QMv l g v KgM v g f v ) 


O ( U ) 


C 9) 


( 10 ) 


( 11 ) 


( 12 ) 


「+ 。（ 1 反 」 3 ). 03 ) 


CM ) 
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粮搪假定，是一致囬的，所以礼定:&，记其特征稂为 

— x ls — h ， ， * ■ ，一九 ” （o < a 】s a 2 < • …< W 则有正交 
变涣矩阵^使 



记\〜 04) 式变为 
八 V 。)） 一 /O ㈣) . 

F (>) — /(，)） 



+ w), 

略去无穷小项，并应用反向 Schwarz 不等式得 


/(^Q - 

/(，)）一 /(#) 〜 




于是我们得出 
定理若 / o ) e A 且有 

— Mgg f < gH(^)g < — r ” gg ’， 

则 ' 

，> — 这 ㈣ _ /( 方 ㈨ ）一 f(x c ^) < (M - mV 

7 ⑻ 一 一 / O ⑻）一 /( P 〕） 、+ ftt /' 


< T 02 • 
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第三章切块法 


述 要 


我们的目的是要在区域 A 中找出函数 

^ Kx ， y ) 

的最大值，在及。中取一点作(仏 y D )， 在（吻，抑）及其邻近两 
点0 。十心 （ q， 少。+々）作试验，得 /(J： D ， 外）， /O。 + A， 
y。)，/U ? yo-J-O 的数值，这里 a , 々以小为好，但又要使诸 
/值可以分辨.作 


(^0 


f(x Q -h h , > 0 ) — y 0 ) 


h 


yo 


yo + O 




K : 。， y。) ) 


过内作垂直于内， h 联线的垂线 4, 这垂线分尺。为两部分, 
丢掉不含有 〆 的部分，留下的部 分用兄 表示，再在其中取一 
点 pi ，… • 

这方法的道理是显然的，因为分割是沿通 过押的 等高线 
的切线，且留下部分包括升高方向的一边.问题在于怎样选 
内？我们有下面几个建议. 

第一个建 议是： 取内为及。的重心，以下将证明，这一方 
法是有保证的.下文将证明，通过重心的任一直线把凸域分 
为两份，其面积的比在4乃及 5 /4之间，因此每次留下的部 
分的面积不大于原面积的5/9，做 n 次后，留 F 的部分不大于 
原靣积乘（5/9)%因而《充分大时，这方法保证是可以无限 

' 10 J - 
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精密的. 

第二个建 议是： 取余下的部分中的法线的中点作为 
下一次试验的内，这方法简单易行，但还没有数学保证，以下 
的附记可供参考.如果力辦出来的实验比办小得多，可以不 
必在 h 附近再做试验，而在内，巧的联线的中点朽做试验. 
如果还小得多，可以在 A ， 内的中点内做试验， 等等 ，得到和 
抑所做出的差不多好或更好笙的点广，再在，用上面的方 
法，这样可以加强这个方法的可靠性. 

第三个方法.假设我们所研究的对象是多边形，其顶点 
各为 G = 1 , 2 , …，/)， 我们可以用 

(+ (^i + , . + + xi)^ +(’【+**. + y f )) 

作为下一个试点.但通过这个点的直线，并没有分图形为两 
部其比例有上下限的性质， 

当推广到 f 个变数时 * 找重心也并不容易，但本章的方法 
并不是没有实用价值的，而理论方面还有深人探讨的必要. 


S 1. — 个几何不等武 ( 二维 > 


定 H 通过一凸区域的重心的任意直线把凸域一分为 
二，其面积之比不小于4/5,不大于 5/4. 

为了替高维作准备，我们在讲几何直观证明的同时，还讲 
易于推广的解析形式. 

Cj 域 D 的定 义是： 如杲两点 （ Ui )，（ Q ， 九）在 D 内， 
则这两点间的直线段也在其中，即它含有 

X ^ t{x 2 — Xt), y = y t 卜彳 — yO 

(0</< 1). ⑴ 

区域 0 的面积等于 


- m * 
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dxdy^ 


D 


其重心的坐标是 


X 


去 I 卜的， 

D 





y ^ \ \ ydxdy m 

A D 



t » 


假定分割的直线就是 t 轴，则重心在1轴上，即 


| j ydxdy = 0, 


D 


( 2 ) 


(3) 


(4) 



(5) 


^轴分£?为两部分，轴上的用表示，轴下的用表示，并 
用|£>|表示 D 的面积，则定理的结论就是 

” CO 


D + 


f l dxdy 

卜 4 

[[dxdy 

D 

5 

11 

5 

JJ 


而定理的假定就是 C 5) 式，也就是 



ydxdy =» — ydxdy ^ 





(7) 
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这 辦是说 如果上半块对 r 轴的力矩等于下半块对 r 轴的力矩 
吋，我们有不等式 CO . 不等式 （ O 的乃' 可颠倒，因而 
得到与定理等价的结论. 

证明分三步进行： 

1) 对称化，也就是在 
保持面积和力矩不变的情 
况下，我们可以把 D 改变 


y 



成为对称于^轴的形式， 
几何直观来 说是： 把每条 
平行 于1 轴的，两端在 D 
的边界上的线段，水平移 
动，使其中点在 y 轴上，这 
样移动后所得出的 W 对 
于 y 轴对称，而且面积与 
力矩都不变. 

解析表示又怎样呃？ 

我们看到，当 y 固定时， x 的积分范围是 


^ i ( y ) < ^ < < P 2( y ) 5 


因此 


\\ dx ^ = \ y >'\ 

= j > s ： 


dx = I 


OfXy) — qpj(y))^y 


i (朽 ( jO — 


dx ， 


以及 


ydxdy ydy j * 






dx ， 


也就是我们的区域 z> 变为 

D’i — (^Pz(y) — 


■ 10$ • 
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< r € —(屮 2 ( y )— 屮 

2 

现在要证明也是凸的 v 也就是如果 （ x 15 %)，（^，>^)属于 
沙，则 

r /( r 2 — a )，y >-= yi +< yz—yO (0 < / < 1) 

也属于后一式是显然的，由于 

一 ^ Cqc 2 (^) — fpj(.Vi)) < ^- Oj(h) 一 ^iCyO)^ 

""y ( 奶（乃）一 <P\(y2)) < h j ( 少 2( 」 v 2 ) — 中 1( 乃））， 

可知 ‘ 

1 

——[(1— 0(中2(力）一奶 Cvi ))+，（ q ^( y 2) — 奶 （ a )〕] 

Ki — OC ^ iCyi ) — < PiC ^ i ) 

2 

+ 心2(少2)—卬1(乃)）]* 

因为 D 是凸的，所以 

(1 — /)<p 2 (y!) + ^i(y 2 ) ^ qt>Xy) 

. (1 — /)^t(yi) -f- fffXyi) ^ qp.Cy), 

由此推出 G ， >0 也属于 O' 、 

因此，并不失去普遍性，我们可以假定 I? 是对 y 轴对称 

的 . 

2) 拉直 .X 轴交 D 的边界于 P ， p 二点，在 y 轴上找一 
点 i? 作等腰三角形使其力矩等于的力矩.我们证明三角 
形 PP 尺的面积不大于公 + 的面积，也就是从 

|| ydxdy ^ jj ydxdy 

推出 
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U dxdy > || dxdy u 


z >+ 


JBk 


三角形的高交的边界于 y ~ y 。， 在三角形内而不在 D + 
的点的 y >>。， 在不在三角形内的点的即得所 

证. 


y 



Q r 


再把和分别延长到 〆 5 使得 〆 7平行于 * 

轴 * 使梯形 PP f Q f Q 的力矩等于的力矩，同样可推得 
梯形 PKp ' j ? 的面积 > D - 的面积， 

因此 

D + > \^ PQR \ 

^ 〆I 梯形 

由于力矩相等， ^ Q f R 的重心仍在^轴上. 

因此，只需证明 

3> 等腰三角形过重心平行底边的直线分三角形为两部 
分 
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jAPg^l _ 4 
I 梯形 5’ 

由于到重心的距禽等于高的2/3,因此三角形 Pp /? 的 
面积是 ^9^ 的面积的4/9,即得所证， 

以下就依这三个步骤，把切块法推广到 i 个变数的情况. 


§ Z /维锥形的体积与重心 


S 


^维锥形的底盘是在 L = 0的超平面上，半径为 r 的 
_ 1维球 

xl H -+ xi-, < r 1 , 

其高为 A , 当 t — 0时，半 
径 P 等于 r ， 当 i — A 时，半 
径 P 等于0,边界的方程是 




P 


h — 

一 k ~~ 


* 


因此，锥形的体积 


v 




dx 


dx^ix\ + * ， ■ 

< r p^x f 




0 


- * dx^xi 

+ “ ■ + < 1 



Jfi 


=3 




!: 0 ) 


dx 


w 




】A 




l h 




这里心-^是 S — 1 维单位球的体积，这公式也可以述为 
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维锥的体积等于底积乘髙的 r 分之一, 
对超平面^ = 0的力矩等于 


卜 ..J 


I* 


h 


Xjdxj * * *dx 




G 


Ml 


I 心 r • •办 


+ 


m m. m 


十 < 〆 -= 似 w ' h 1 \ o - 0 ^\dt 

0 


= w^r c ~ % h 2 is^s -H 1) + 
因此，重心的坐标是 


^2 




0， ^ 


h 


■ 


s 


过重心平行于 I = 0 的超平 面把/ 维锥分力两部分，其 


是高为原高的 


f + 


穩 ， 即依顶点缩小 


S 


S 


的 


圆锥，它是原来圆锥体积的 


S 


S 


倍，因此踽锥被过重心的平面^ 
积之比是 


h 

■^圓 

$ + 


分为两部分，其体 


S 


(-s-i 


这相当于上节的第三步 




S 3 .对称化 


我们假定 < 维凸体 D 被过重心的超平面 

0 

分为两部分，上、下半邰各以£)% ZT 表之，则问题变为在 


m 


_ 
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x s dxi * * ■ dx 备 


- j …卜, 


dx t " " h dx s 


u 




的条件下*求比例 


* _ ， 


u 


■ 

dx^ - *dx s ^\^ * * 

" D~ L 


i/jf L . . ' dx t 


的上下界的问题* 

先用以称化的方法， 不妨— 

假定乃是旋转体，再由“拉直” 

法来得出我们的结论 . 

为叙述简单，我们假定< = 

3,将 O ，， a , A ) 写为 O , y ， 

■ ■ 

依 w 平面对称化，一垂直 
于^平面的直线交凸体于二 
点 p ， L 把这线段的中点移到 
^平面上，则得出另一体认， 
D ^ Dt ’ DT 的体积和力矩仍然相等，即 


S 


團50 


11 ( zdzdxdy ^ | j | zdzdxdy y 

〆 . D i 

j zdzdxdy ^ jjj zdzdxd ? ， 


D 


D 


dzdxdy Jjj dzdxdy ， 


d 十 


D 


J J I dzdxdy = J I ^ dzdxdy m 


D 


⑴ 


V. 


( 2 ) 



⑴ 


这一点的证明是很容易的，因为当定了，7的积分范围是 


_ 


m 
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cpOr ，^ < <p(x ， z) ， 


HI fC^)dzdxdy j(js)d^dx 


« Ka ，） 

^ CsTtf } 


dv 


因而可用类似于纟 2 的方法证得 （3) 式，并且同样方法证得 
0 X 也是凸的， 

实际上 3 以上的方法可以用任何一个过=轴的平面来代 
替，这样我们得出一系列的凸域，这些凸域的极限是一个绕= 


轴的旋转体，为了证明这一点，我们讲得更具体些, 



图51 图 52 

先固定一个 a 截凸体得一平面凸域的周 
界距原点最远、最近的距离分別为 C 与 C ， 点在 P 与从 
原点作直线垂直于，， P 的联线，依这条直线对称化，得出来 
的图形长径一定比 C 小，短径一定比 C ' 大(注意，由于等腰三 
角形顶点到底边上一点的联线一定比腰长短，到庇边延线上 
一点的联线一定比腰 长长， 因此对称化总是把 C 缩小，把 C ' 
增大)_ 

由于得出一联串的 C 是递降的数列，因而有一极限 
一眹串的 C 是递增的数列，因而有一极限 C :， 如果户 C ;， 
则以上的方法还是可以做下去,因此所得的极限是一个圆. 

已经成为圆的对称化也还是它自己，再桧查任意 I 这做 

■ ■ ' 
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法可使任意 Z “的截体都是圆. 

但需要声明一下，一系列凸体的极限还是凸体，因而证明 
了：存在着一个绕=轴的旋转体，适合 C 3 ) 式. 

S 4* Brun-Minkowski 不等式 

对一凸旋转体 来说： 怎 〜所截得的圆的半径是 
面拐是 ^ p 2 (^) = /(«)，因此 

是一个凸函数，也就是有不等式 

(l-a)p(a) 如办 0 ) ^p((l-d)a-h3^)^ OOCl w 

对£维来说，^ = a 所截下来的 Y — 1 维体积以 V(a) 表之， 
则有不等式 

(1 — -h 3V{h)^ 

< f((i — d)a + 

(0 < & < 1), 

这是有名的 Erun-Minkowski 不等式 C 见 Minkowski Geometric 
dcr Zahlcn , 数的几何）. 

另一个 办法： 先证明了这个不等式，然后以*为中心，作 
一 体积为 VU) 的 i 一 1 维球.当〃依 2 轴移动，则得出一 
个凸旋转体，而凸旋转体适合于上节公式 （ 3 ) 所示的关系. 

S 5 .拉 直 


主要定理 通过/ 维凸体 D 的重心的超平面把凸体分为 
两部分，其体积之比为 A /， 
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事到 如今， 这个定理的证明就不难了.由前已笳，只要假 


定这个 f 维凸体是重心在平行于底的超平面^ = 0上的绕 
t 轴的旋转体 ，作 一个以1=0交0的| — 1维球为底的^ 
维锥体，它的高 A 使这个锥体的力矩等于 ZJ + 的力矩.设这 
个锥体交 D 于I = G 由于外锥内的点适合于 o ~而锥 
外内的点适合于 x s < h 因而锥体的体积不大于的 
体积，把锥延展到下半空间，使“台”体的力矩等于 D + 的力 
矩，同法证明“台”体的体积不小于 ZT 的体积，因而 


m 体积锥体积 

D - 体积/ “台”体积 

■ ■ 


即得所证. 



5 &说 明 


1) 这个方法有很大的优越性，每次留下的部分不大于上 


次的 


倍，由于 



所以留下的不大于上次的 
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倍，这个数与 J 无关，而每次要做 f + 1个试验，因而如果要 
求使留下的体积是原来的那末试验次数的数置级是 

m 

jlog 

即使用体积是长度的 f 方计算，也只要 

s 2 \ogtn 

数量级的次数计算就够了，这也就是一维的？倍. 

2) 这方法有两个缺点 u —个是差分 

IM m I - ^ ■ I 1 ■— m 

h 

代替微商，“小数除小数' 这是计算数学中忌讳的，必须小心 
在意 ■ 

另一方面除^ 2外，重心的计算很复杂.关于这一点 
可以求助于电子计算机，通过求积分的近代方法，问题可以化 
为球面上求积分来解决. _ 

n gL . r n v 忧 r > 山片一占 p _付 汶占柞 1 ^ 一 
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第四章二次通归法评介 

述 要 

以下介 绍一个较简单，较粗糙的方法，虽然如此，但比一 
般利用统计学上二次迴归求最大値的方法还可能好些. 

方 法是： 作两条平行线 A ， L 在这两条线上用单因素优 
选法各找出一取最大值的点 Pt , P 2, 通过这两点作一直线济, 
在 m 上取最大值的点就作为优选点了，所取的值就作为 
我们所要找的最大值. 



这方法的局限性很大，但还是有些启发. 


S 1 •背 景 

正是从椭球 的等高 线图形-簇以制髙点为中心的同 

心椭圆启发出了这个方法.它的基 础是： 两个同心椭圆的两 
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条平行切线的切点联线一定通过中心. 

其证 明是： 假定中心在原点，两个椭圆各为 


这里 


ax 2 4" 2bxy + cy 2 
ax 1 Hr 2i?xy + cy 1 


彡 2 一 0 . 


=； 








他们在 O u )0, 以的切线各为 

-h cyj)y ==* l 3 

(_ax % 4- byi)x + (Sx 2 -h cy 2 )y — 

如果他们平行，则 


即 


aXj "h bvi - hx^ + c m y x 

ax 2 hy 2 i^x 2 + cy 2 ^ 


— V)Oiy2 — ; 0. 

由于 w — 炉尹 0 5 所以 


£j^ _ £2 

yi 


也就是 （〜 yj ，0 2 , y 2 )，（ o , o ) 在一直线上. 

本章开始所讲的方法，实质 上是： G ，込是 两个依中心放 
大的椭圆的切线的切点（并不排斥两个椭圆重合的情况），因 
而他彳门的联线通过中心.中心既是这联线上取最大值的点， 
也同时就是制高点了. 

这些结果还可适当推广到稍为广一些的函数类.假定 

<0 fO ， y)<l 

是一凸域，并且适合 

〈 2 〉 yt) = y). 

上面例中的同心椭圆族属于这类函数，又例如 

f(x y y) = V k| rt + \y\ a (a > 1). 


* ii ? ， 
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满足上述条件 <1> 及 的 z = / O , y ) a 等髙线由 

y ) = 1 

所形成，这是一批以原点为中心的凸曲线.从中心出发的某 
一方向 

丨 =/ cosdy y = / s；n 6 

上的一点/ = k 处，曲线 

( A 0 = ^/( cos 0, sin By ) 的切线斜率是 

一 i / 空 L — + 邕逆 V 

dx / dy \dt dx dd dx // 



$ 十 H ) = - [以 0 sw>os & sia 


sm 


9 /(cos ^5 sin6)] ^/[p^\n 0/(cos0^ sin9) 


+ cos B {C cos 0 5 stn ^)1. 

50 

这是与 > 无关的，也就是在这一方向的任何一点上等高线的 
切线都是平行的，这是上面讨论的一个推广. 


S 2.二次迴归 

这个方法的限制性是极大的，但与统计学上二次迥归的 
方法比较，并不见得差.下面介绍一下统计学上的二次迴归. 

我们预先假定所考虑的间题合乎二次模型，也就是假定 
、是 r ， y 的二次式”，即 

z — ■ = A(x ^ r*) 2 + 2fi(jc 1 ― 1 x*)(y — h 广） 

+ c(n ， (i) 

'，，，，，人 AC 是六个待定常数，如果我们在六个点上 

* m • 
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作了试验，得出六组％，= 1，2,…，6)，六个待定 
常数就可以确定. 

统计学者把试验中可能有误差的因素也考虑进去，他们 
试验爪次，在处试验出来，然后求 
A , B 3 € f 使误差平方和 

W 

3> - S ㈦ ，一 — 2B ^ 

i - 1 

一 ^*)(yi — y*) ~ — y^) 2 Y 

最小. 

由 （1) 可知，如果 J < 0,且 ac < s 2 , 则显而易见在 
T = ，处，=取最大值？ . 

这个方法的计箅比较麻烦，而且一开始就来了一个限制 
很强的数学模型 （1) 这一假定， 

这模型的等髙线是以 W ) 为相似中心的椭圖族，反 
之，等高线有相似椭圆族性质的曲面不一定是 Ci ) 的承式，因 
而，本章所研究的方法虽然特殊,但比二次迴归还较一般些. 

预定数学模型是二次的，这的确是一个极为大胆的设想, 
它很少可能符合客观事实的，但是当接近最优解时，也就是 
^ 7 ，充分小时，我们略去髙次项，不妨就用 n ' 

y 的二次式来代替 L 而这正是难以辨别好坏的时候. 

§3. /个因素的问题 


先讨沧三个因素的 情形： 

开始取二个平行平面，在其中一个平 面〜用 上面的方法， 
求三条直线上的最大值，得出在这平面上取最大值的点.在 
另- 个平面 A 就不需做三条 S 线上的优选法了，因为前一平 
面已定出了二个共轭方向 A 与纟3.我们在平面巧上取一条平 
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行宁[的直线 C 在其上优选的最好点过作平 行于％ 
的直线/;，在其上优选得最好点则 P ; 便是 平面〜 上的最 

好点了,再在连结 与〜 
上最 好点^ 的直线上优 
选，便得出所求的最好点 

因而三个因素时，总 



共做了 3 + 2 


6批 


单因素的优选试验，然后 
得出最大值. 

我们很容易用归納法 
证明， y 个因素需要傲 


0 


2 


批单因素的试验. 

如果单因素的试验用黄金分割法去做，次数的数量级仍 


然是1叫丄，现在只是前面加上一个 


1) 


2 


的 因子. 


§4.讨 


虽然屡次指出了这方法的局哏性，但这个方法本身还是 
有所启发的+ 

首先，因子7 — 1实在太大，我们可以考虑以下的 方法: 
从一点 y t ) 出发，联系切块法的思想，作一最陡上升 
的方向 

L ^- 丄 : M， >、) 一 f( x ， …) tt 

\ h 
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h + Ka ， : v! + 走）一 /Oi ， yjy 

垂直这个方向作一宣线卜在 G 上行单因素优选，得一点 o 23 
h )， 再在 O ,， y 2 ) 的联线上找最优点 i ?. 



这方法推广到£个因素， H 要 y 个单因素优选法就行了 # 
其次，我们丼不拿作为最后数据，而作为再一次的开 
始，闬这办法来逐步逼近最优解. 

其三，改变迴旋上 升法. 

在一直线 G 上找到一 最优点 
Pi ， 通过&作垂线 I 在这垂线上 
用单因素法找到另一点过^ 

作 G 的垂线 L 在： 3 上又找到点 
&，接下去，不作 匕的 垂线了，囱 
于 P ” 5是两等高线的平行切线 
的切点，因而联 p 15 / V 在这直线』 4 _ % 

上找最优点/ V 过5作“的垂线 A 等等， 

从多 因素优选法可以看出，随 着维数 r 的 增加试 验次数 
很快地增加，所以我们切不要因萦愈多愈好而轻易地增加变 
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量, 而应当抓住主要矛盾 ，先 集中少 数几个 因素来做试验.这 
种突出主要矛盾的方法可以大大减少工作量.当主要矛盾解 
决之后，另一个矛盾上升为主要矛盾时，我们再用同样的方法 
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第五章抛物体法 


述 要 

单因素抛物线法是利用三点的数据作出一抛物线，以达 
到此抛物线最大値的点作为新的试点，这方法所需次数的数 
量级是 


1 叫 1 叩？ 

把抛物线法的想法推广到多因素，便得出抛物体法. 

以二个因素为例，拫据在〜，卜1，2 5 ■■*， 6) 
点上的试验结果，可以确定一抛物面 , 

z ^ y) = ax 2 -jr l^xy + cy 2 -h ^ -h + f # 

求出此拋物面的顶点 O % y % 然后在 O % 作试 
验，与&，匕 f %， 心匕的结果一起又可构成新的抛物面 5 …. 
这样不断叠代. 

我们将征明，当满足一定的条件时，这方法收敛速度是相 
当快的.在/个因素的情形下所需次数的数量级仅为 


s 2 log log — , 

e 

s i . 矩阵符号 

为了便于讲述本章的方法，我们引用矩阵符号，以」= 
表/行 m 列的矩陴； 
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a == (^) (i < < y < m ). ( i ) 

向量 A ： — a , …，， &) 可以看成为一 行《 列的矩阵 
V 1 '我们常用/表示由^经过行列互换听得的矩阵，它是 

馆行/列的. 

定义我们称 


Mil - S S f^I ⑺ 

1< I <1 !<?<*» 


为矩阵 W 的模.特別对于一个向量X，我们用 〗x!f 表示它的 


各个分量絶对值 之和. 这样定义的模显然有以下的性质: 


\\A-hB\\<[\A\\-h\\B\h (3) 

14 B \\< U \\^\} B \1 (4) 

当然，这些式子，只当 A -^ B , AB 有定义时才有意义. 

引理1 命 X = X U _U = 心，若 

1, ⑴ 

则 

||X(/-^)-1 < l/(l-t,)|lX|f ( (6) 

证明我们有 

J 1 ^ (/ — A 1 — ^ C.I + 3 + * - * )， 

由 （3) 及 （4) 可知 

\\X0 - jr x \\^\\X\\ + \\XA\\ + \\XJ 2 \\. + •■■ 


< || x||(i + « + ¥+•••) = l/o — WM , 


S 2 .方法的背景 

我们的问题 是：求 j 个变数 r 仏••*，&)的函 

数中00的最大值. 

在丄（』+ 1)0+2) 个点 
2 
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，)， ⑴ 
V 0 〉+ m (1 (i (s )， 

方⑷ + ”(。 + q) ( 1 ^ ^ ^ i - ) 

上各做一次试验 3 此处 

q = (0 ， … ， 0, i ， 0 ，•，，()）， 

、 ■■垂 •— 〆 

f 一 1 个 

得出丄 0+ 1)(，+ 2) 个 数据： 

2 

少 ( ：⑻），子 0 〔 0 ) 十的）， 9>(>(° } + ”（~ + 〜））• （ 2) 

我们现在找一个二次式 

少（戈 ）= a + 2bx r -h xCx\ (3) 

其中 《 = b - b ^\ C^c ^ c ( ^> 为一定负矩阵，使 
在 （1) 中各点上有 

400 0 (fix') ⑷ 

成立，现在总共有+ 0+ 1)0 + 2) 个方程，由此可以定出 

▲ 

4 (，+ 1)0 + 2 ) 个未知量 〜 h C . 

2 

令 

A//U) = -i- [f(x 4- ^ei) — f ⑺]， (5) 

则 

A f( p(x (o) ) - = — [2 办 (> ⑻ + 的 y — it{x^y 


及 


+ ( 尤 〈 0 > + + ”<_■)’ —- x^Cx^] 

^ 2bc\ + Zx^C^t 4 - ” qCe; ， (6) 

A i A , q &( jc (0 O = 2 e f Ce ) (7) 
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由（7〕得出 


由 （ O 得出 


( A ，( pO ⑷)，…， A 冲 O ⑼ )） = 2 ^ 4 - W 

+ | 刀 (AftpO ⑷) ，…， 你(，))） 

即得 

■ 

2 厶=! — — Ap ( p(^ (o> ) J * * * J 

•2 

L 

△冲 (尤 ⑷） 一 + 巧 △>(> 忉 )）一 2 t Co ) C + ( 9 ) 

-U 


再由于 C 是定负的，及 

少 O ) ^ <X — bc_ l b’ + O + i ? C ^) C{x 4- be - 1 y 

所以当 5 


x _ 方⑴ = 一办 c 一 1 

时， 0 W 取最大值，因此 

太⑴=/ 0 〕 — t ， 0 〈°))M (>〕）-、 

这里 


7( 方⑼） = (ji ( cp(^ 0 ^ — Jl_ 7jA^p( ； r(°)) ， * * ■ ， 

■ 2 * 


( 10 ) 


■ ApO ⑻）一丄 V A ]< p ( x ^)), 

M ( ’卜（厶為 W %。, /<,• 

我们的方法是利用递归公式 

方 ㈣ ） = 士 ) — V ( x < ， r )M ( x ：，yyi ⑴） 

来逐步逼近 vOO 的最大值，此处 

^(^ v} ) (△』 〆 #) 一丄 ”△>(#)，…， 

A 史 O (’0 — + 々△->(>“〕））， (12) 
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及 


M -= ( A ; △冲 (>“)））!《, 如 (13) 

假定这最大值在 a : = / 处取 ，则 

■^< pO ) = 0 (，■ = 1， 2, …， j ). (14〕 

为简单起见，以后把 I ^ Or ) 记作 9 ^. 

在证明这力法的有效性时，我们要用到下节两条定理. 

■ ■ ' 

S 3. 两条定理 

定理1如果 

S s 

d SI < P^I < 夂， (15) 

r = t i=l 

则 

一 /I < 丄 K Q ]^- i } — 太(， 

+ 2 Vl \ x ^ - / 1 )|| M (汐 m (16) 

证由 （ 11 ) 式得 

( V 叫一 = — V ( x ^) 3 

G ⑷一 = — FC〆 1 》)， 

因此得 

0 (1J+1 ) — x ( v ： ') M ( x w ) V ( x iv ~ iy ) — F (. r ^) 

+ ( x <v) — AKA -”)， (17) 

注意到 

^0 dtx 

逆 Xj ^ x< ^ } — i ^^( 尤⑺ + 0P(VW — x iv} ))da ? (18 〕 

%=1 
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t 


及 I 

△iK —、） 一 fxi{ xiv) ^ 0^/7/Vos 

Jfl 

(17) 式右边矢量的第；个分量是 

(△/ - 合 厶】 )(<pO u_1 〕 )- 屮 0(*°)) 

+ i ] ( 4 V) - 4 V _ ”) 厶，△沖 C - 〉）• 

-S 0” - 々)){(△‘-★ ㈣ ) j : T ^ iv) 


+ c < W -” 一 X^)}da 
一 L 9^(々一) + (XTj€ k )d<i^ 

-X wr ”— 中 〕 { I : j : 

A ^ 1 

4 - ot(;r (v 舞 ” 一 々〉） 十 ^Ji€i)dad^ 

—+ ” L L L 叫龙 /〆 10 十 a(x (v ~ l) — x iv) ) 


+ + y ^ e ^ dad^dy 


s: s: 


x 


Cn— o 


2r\e^ H- P7]e：}^adfi 




S (^- l) - ^ v) ) 


k 


m\ 


0 


[(1 — a) ( 泎 ) 


4 一 )) - 

£367 /4 + pi]e i + r((l — 0 f)(x (1 ° — A 一”） 


ctije ^] dad^dr 


2 


V 


m >( 


X 


(w 


-+- a(JC (y —) —— 4- Ji]Ci)dckd^dy | 


■ 


u $ 




(19) 
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其中 




由假定 （ 1 5 )， 可知 


一 H ，叫卜 f ^|)^ 


iO 


X 


(: y —1) 




fl 


.{T 「 (1 

Jo jo Jd 


^ dad^dr 


Ji 


0 ^ 


adad^dy 


j~v& 。 dadpdr j ||M(jt c ^) -1 !| 


2 


( If ^ 


!•■—1) 


X 


< y)j 


2 


2*(-1) 一 m!A/W”|- 

这就是 （16) 式. 

定理 2 如果 

4 




及 

则 




r 


< 1, 


IlMCr^tl ^ 


1 


a 


|| M( yH 


V 


证命 


则 


M (’>) 一 M ( x ^ y )= (⑹）， 


Pii " △^△ ，（史（ 〆 10 ) — qo(r 卜 ”)) 


厶 _△: 


o da 


<p [ r Cl3 


-ij 


+ a( x 


( 20 ) 


( 21 ) 




x 


Cv 


ly )Ucc 




t OP — 十°)八為 


fr - 


l 


- \ q>, k (x u - l) + a (x 


f ^> 
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因此 


J / ^ 


\u - S 2 i^i * i^i 


i^=i 


<t±±±\ 

^ - t / = 1 1 f 




x 


r ij i 


# 


o 


△ £ △沖々 卜 !) 


ae(x 


C ^) 




命 


一 / 卜 i)|j - f|M(^‘ 1 ))-】||. 
/ — M 


學 


则 

= M(x^- l >y\i - Ay\ 
由假定 (20) 及引理 i 得出本定理. 


( 22 ) 

(23) 


M •有效性 


则 


定理3 如果 

欠 & ⑴一 


( = 0.382 …）， 



j|^ F ) „ ^-i)[f ^ 2jh 


\ W ^- x ^\\ 


< 



证 


0 由定理1，2及 （25) 可知 
K \\ x ^ - ^ ( ^||.] m (^)^ 1 !| 

< (K\\x^ - ^- l >i[||M(^)^ll) 3 



(24) 

(25) 


(26) 


，東 30 • 




h-^c 
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现在用归纳法证明 




0 < 〜< 


3 — V "? 


2 




(27) 


由于 /0 O 


X 


2 



X 


x 的递增函数 (0 <^< ]) ? 所以对 


于不大于 


/ U ) < / ( 


/3 


V 3 


的#可有 


3 — 


2 


h 

)- 


3 



5 


2 


/( - 


+ VT 


2 



2 


V ， 


2 



V ’ 5 


2 


— J 飞 

~ 




于是 (27) 得证， 

2) 由定理2,可知 


a 




V 


< 



2 


i ) 


< 


\/y + 

2 




IlMWjj, 


(28) 


再用归納法证明 （ 26 ) 式•当 V = 0,此式显然成立，由 
定理1， _ 

I 卜 ㈣)— H < 到汐)一 w- 1 ) ㈣ AJ0 C ， ] 


<K 


3 



5 


2 






5 — n 


2 


琴 


IM 1 ) 




i ." 


< 


3 + VT 


2^-2 


2 


«D 



2 


5—1 




2v_ 
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• Ik ⑴- ， il ( 〆 5 - , 」)、. 



即得 C 26) 式. 

定理 3 说明，一旦找到一个初姶值适合于条件 （24) 之 
后，我们这个方法进行得是非常迅速的，其收敛因子是 


第一个是主要的，第二个因子和黄金分割法差不多好，而且 
维数增加，试验次数只增加 P 倍，而不是 f 次方，如果联 



时，就可以使 

||，+” 一 汐>|| € ⑴一 ，I 

注意：在 

| 卜00 一 X (V - 1：> \] < 27} 

时，我们停止做下去，此时，我们已经尽可能地达到，的渐近 
值了， 

附记也许有人提出条件 CM) 不易验证.我们的建议先 
用其他方法做到相当程度后，再用此法，在用此法时，如果试 
验出的？>0巧比 中 (#) 大，则照样做下去.当 < p (^ ay ) 不比 
K，） 大时，可以与切玦法同时考虑，因为在，⑼+从, 
^ --H nOi + 做了实验，可以来检查， 1) 
个点)的上升方向，而另外0 + 1) 个点，+畎又有 
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了， + 1 个新方向 j 艮据这些线索可以得出棱小的区域，但有 
一点理论上是肯定的，充分多次之后，这方法便一帆风顺了. 

S 又补充方法 


在多变数时，以下的简单方法也能起到好的作用.我们 
用叠代公弍 , 

(1) 文 ㈣= x^-u(x^)M(x^y x , O = 1 ， 2,…) 

也就是 Af 并不变化， VO ) 代以吏简单的 

W =(叫，叙2，•，‘，4)， 

— 〜 0O 冲 ^qo(y), (/ « I ， 2, “ 、 J). 

类似于前述，容易证明这时有 

\\ x ^\- x ^\\ < K \\ x ^ — ^-^[| (| k Cv ) — 

+ 丄[| 2 ⑷— 21 + 斗 \\ M ( x a y % 

2 2 / 

为了便于比较，我们略去了 ^又记^为制高点， 

■ ■■国 ■ ■ 

概括地说，抛物体法每次要做 

丄 （J + 1 )(V + 2) 

2 

个实验來定出 

中 O 。）， 免 Oa 十 O, <p ( 尤 0+ V “ + . 机 X 1 < ^ < ； < ^) 

的数据，而收敛情形是二次的，也不难证明 

| 卜 ㈣ — 《 || 少 ) 一， |f 、 

这就是说，收敛指数成倍地增长着. 

补充方法是每次做 f 十1个实验来定出 

屮(尤）， tp(^x -h 7]“) ( 1 ^ f ^ /) 

的数据，但收敛情形是 ^ 

Ik ㈣ 一？ II《 P 
即收敛指数上升一阶. 
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第六章与计算数学的关系 


§1. 问题的叙述 


如果 / O , o 是一个有了表达式的函数，则求 

a = fO, y) 

的最大值的计算问题，和我们所讨论的问题一样，但"试验”巳 
经可以更具体地理解为“给了 t y 的数值，通过计算得出 Z 
的数值”.由于近代电子计算机速度很快，因而多做呰“试验” 
是完全可能的.如果计算一个7的数值需要一千次 3 在每秒十 
万次的计算机上只不过1/100秒的时间而已.但如果我们的 
h Y 各要求算四位，且对每一对（心 y ) 箅一心，需运算千 
次，然后比较出最大值，这样总共的计算次数就是 

1 0 4 X 10 4 X 1 0 3 10 1 、 

J m 

即一千亿次，这样每秒十万次的电子计算机还要计算三百个 
工作 A 时，太多了： 

如果变数多了，增长率更快，电子汁算机虽快，但还是不 
能胜任. 

工作可考虑分两方面来进行，一方面搞电子计算机的往 
加快速度方面努力，使每秒钟的运算次数愈多愈好，愈可靠愈 
好，另一方面搞计箅的，从减少计算次数方面努力，想方设法 
来减少计算次数. 

前面各章所介绍的优选法，实质上也就是求最大值(或最 
小值)的快速计算法.实际上不仅如此，很多计算间题，都可 
以化为极值间题来解决，而极值问题离散化后，就与本书所讲 
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的方法有关了，因而优选法在计算数学上也有其重要性. 

例解超越联立方程 

Y) — 0，《0, y) = 0， 

可以化为求函数 

y) m f + g 1 

的最小值的问题， 

所以，我们不能粑优选法简单地理解为一个减少试验次 
数的方法.实际上，这个方法给计算数学也增添了内容. 

& 2 . 黄金分割法的计算格式 


求方程， 

y /( a ?) = Q 

的解.假定在幻中仅有一解(或本法仅能找出其诸解之 
一 ）4 十算步骤如下： 

1) 写下 b , w 5 - 1 - = Q .6180339887; 

2 

2) 编出计算1/00丨的程序； 

3) 箅出 

a w(J> ™^ q b —— w(^b ——心 l 

4) 依 2) 算出 |/( p)U i / WI ； 

5) 如果在合乎要求的精密度下1/001=0,则停机并打 
出 P 的数值，同样处理 IK ^) I ； 

6) 如果不合于 5) 的情形，则比较 I / O ) 丨与 I 關的 
大小：？如果 lf ( p )! ^ I f (今 ） I ，则取^ M p ; 如果 

Ifip) 1 < 1/( 彳)“ 则取 h 工 q ， b\ 

7) 如果么 一 q 合乎要求，则停机打出 〜匕及 I fM \, 

IKMI ； 
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8) 不然，回到原来的程序，但以 匕 代替~夂 
注意： 如果在 5) 停机 5 我们找到了近似解，如果在7〕 
停机，而 \ fM \ 不是合乎要求地小，我们可以得出结 论：不 
是问题在 （ A ， h ) 中无解，就是机器出了毛病，从 [与 
i /(^) i 的大小，可以来纠正错误. 

如果 /0 O = 0不止一个解,但能估计出每个解所在的范 
围，用上法可找出所有的解. 

§3•对开法 


对开法的程序也不复杂:我们可以用来解联立方程 

y) = o, y) -= o. 

假定在 ^ 

n 4 x 《 b, c ^ y ^ d 

的范围内有一个解. 

i ) 用黄金分割法算出 

F (*，含心 + a ) — f ( x ，+ ( c + <)) 


+ g ( r ， 士 O + d ))!， (；2 < Jr < b) 

的最小值 h 在 处取这值，如果在合乎要求的精密度 
下，则停机打出 A ， + (c + 

2) 同法算出在 （ + 0 + 6)， yi ) 处 P 
取最小值^如果〃也合乎精密度且7=0,则停机打出 



3) 分四神情肜 
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当尸 > ‘7, h > 士 （ c + W 时，取 

2 

6 士之〜 = 丄 O 十 d )，< = A 

■dW 

当户 >《， yi <~( c-h d ) 时，取 

a l = ^1 = C! — i' ? ^ = A ( c + J); 

2 

当 p < l 々>!(> + f ) 时，取 

•2 

a i = ~(^ + by = c x -= C, d x — d\ 

当 p < 《， a < 丄 O + 6) 时，取 

〜 + (tf + &)，G ™ r J ^! ==» d % 

4) 以〜， 1〜 A 代 替〜匕 c，< 重新开 始. 

附记 h = 时’ P 不能大于？，因为在7 = 

+ 0 +心上不可能取比 f 更小的数值 t 因而本 法仅在 

p U =音 （ tf +从 yi _ + (c +乃对失效，也就是无 

i 仑横向或级向， F( Xy y ) 都在 (+(« + 0，+ 0十心)处 

取最小值，为了避免这一情况，我们在 2 )后 U 前加上 一条： 

如果 f =亨， + + *)， yi _ 丄心 十心停机打 

^ 2 

出^ g 之$等数值， 

这点可能就是解，因为4 0与 /( ) _ 0 

dx dy /v ^ y/ ’ 
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^ 0 等价. 


5 4 .旋升法 


与 5 3 同样的问题，旋升法的计算程序 如下： 

1) 用一维方法算出 F ( x , 0在 

a y = yi ^ c 

上的最小值，并设在 r 1 时取此值.如果此值合乎要求， 
停机，打出 A ， 

2) 仍用一維的方法算出 ）） 在 




c ^ y ^ d 


上的最 小值. 设在 y p 时取此值，如果此值合乎要求，停 

机 3 打出^ y 2 ; 



3) 以竹代替 h 
重复以上的计算 . 

这一方法实际上和 
以下的方法同一 原理: 
先录 /( r ? yi ) = 0的解 
^ =" ^1,再求 >)= 
0的解 y —外等等.请 
读者注意，把这方法用 
到二直线 


(\ p \ > kl ) 上去的情形 • 若先求 h = F 的解^后求> = 
px , 的解 y I5 •.…-，便将偏离二直线交点愈来愈远.但若兜求 

的解〜后求 y ^的解1， . ，则能收敛于二 

直线的交点.从这一个简单的例子，可以启发出一般原理来， 
甚至可以启发出收敛的条件来* 
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S 5 .数值微分法 

■ ■ ■ ■ 

■ ■ ■ ■ 

我们在研究单变数时 5 往往用一点及其邻近点的差分来 
逼近微分，也即 

补 一 & 

在多变数时，必须增加另外的假定，而不能泛泛地用一点及 
其祁近 几点的数值来逼近偏微商，例如在$ = r。，y D yo 处 

y ) 的偏微商 H 的值，能否由此点及其两邻点 

Ox Oy 

(^，”)，（^， yi ) 的函数值来逼近呢？不增加条是不行的, 
譬如说，三点在一条直线上不行。因为这样只能决定一个方 
向; 夾角太小也不行.另一方面，如果跗方向趋于 Oo , 的 
速度太不一致也不行.例如 

/0, >) = JT 1 + y 2 

在 X = d， V =« 0时偏微商都等 f 0 ; 但^ = £，卩= 4时 

-4 - \}^ £ 文 -4- g 1 

lim -- 雄 litn - -= 00 

x * e 
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图 61 

这是飄第五) 章中我们不是随意取…点 (或 > +1) 


0+ 2) 点）的道理， 

对怎样的三点有此可能呢？实质上，如果面积的絶对值 

■ 

t …门。卜 (彳(-…二一# 

l x i 一 h yi — 


+ \/(n — 丈 。) 2 + (n — y^) 2 } 2 . 

在此条件下，可望成功，读者自证这条件保证了夹角不能任 
意小. 

为了方便起见，还是采用矩阵符号 ，在/ 维空间给了 r +1 
个点 x u> 0 D 0，1， ■ * • ， o 从 

<p(〆*)) ■= a + bx a) (i ™ 0 5 1, 2 5 (1) 

定出向量6,在什么条件下，向量 6 可以逼近 



定理 1命4 
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如果 M f M 的最小特征根 > 欠乏 ] d)M 

o — gradqpO 〜 ))0 — gradqoO ⑻ ))'<,Q ^ d), 

i - 1 

这里 0 仅依赖于 9 的二阶微商的上界及人 ' 

证由于 

<pO ⑴） 一 甲（ : r ⑻)= f -^<pC^>4 - a(> ⑺一 x^))da 

r Jo da 

= 2 (//> — xf^) T <p Xj (x^ + a(x (i > - x^))da 

J -l 』 0 

J jr 

= 2 (4° — 4〕 > )(^(>(。)）+ D C4 ；) — x f } ) 

i-l / =1 

- \[ — W))) _ 〜（，）] 办 

0=5 grad 史〔欠⑼ ）（ 尤 “ ）一 x^Y 

+ 22 - 彳))（中 — 4°)) 难， 

; k 

此处仅依赖 T P 的二阶微商， 

由 （ 1 ) 可妇 

⑴ )— （ ^( 太 ⑻)=— : 

囚此 

(b — ㈣ <po co o)o ⑴一 = p t 

这儿 J 

内 = S S w° — -rx^ y - 逆)也 

•i fe 

及 2 

厂一常数，」仅依赖于 p 的二阶微商的上界 . 

* Hi 


(3) 
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( 3 ) 式可以写成为 

O - grad - ( ?15 灼 ，…， 化）， 

因此 

O — grad 史 O f0 )))M’M(> — grid tp(x^)Y ^ ^ ^ 

由定理的假定可知 

K 名 , i(J> - gr 今 {> ⑼ ))0 — gr^d<f(x^>y < d) 

此即所求. ' =l 


§«• 方程组的数值解 


本书第五章所介绍的方法也可以用来处理來方程组的数 
值解的问题.求 




A 


的数値解，用矩阵写法 


) M h -，，， L(〜xj 


0 (1) 


fM 




m 


先给了精密度7?，例如精确到五位小数，则眾习 
用叠代法 


X 


t 儿 


c^+n = 


iAxJ(x<^) ^ 


( 2 ) 

10 —， 

( 3 ) 


A(x (v >) 


△W") 


ai/oo = ./u + 


_ 


( 4 ) 


定理 1 假定 


n 






dxfdx^ 


< 


( 5 ) 
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则 


一 < 丄 M ( I 卜卜” — 

2 


+ dip -” 一 

证由 （ 3 ) 可知 

0 Cs+n — x (v 0 <^ Cv) ) - 0(^- 1} ) - / O ⑷) 

+ ( x ^ — x ^-^ ACx ^) 


—\ y- f( x<：v} H - «C 

Jo 


x 


c ^- l ) _ 


x uy ) d & 


fi 


2 ( W ) — 


fj 


1 = 


2 os 


w — l ) 


々)） 


0 


Xt { x ( ^ f - 心卜” 


i 


#> Va — I 。 f Xt (x 


Cp -1) 


77 


Sc 


X 


f 卜 1) 


中） 


0 


1 d 


j d 浮 


aljjj)da 


，心卜 ” + al 


iV 




K(l — c()O <t0 — / 卜 ”）— cc^r})]^^a 

j 1 (^[Cl - a) 

i -1 j0 jD ft ; t 

(;— o — ⑹划飞 m〆 咖 


* 


因此得到 


n n 


I 


汐， +1) 一汐) 


<{ a / I ] Sin p> i 卜 


IP 


Cv— 


X 


l) ! f [ (1 一 ct)dad^ 4- tjM U |jt(,) 
J0 J9 / =1 

°l \[ 叫 j 


K ^-(\\ x ^ 

2 
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# - if + y 卜⑻— 


定理 2 假定 


n 




d 2 h 


OiTjdx^ 


< M ， 


M 




⑺ 

(8) 


则得 


M Wll € 




[If ”）1 


( 9 ) 


证命 

A{x iv ^ - A(x^)^ O") ， A{x^) 




(如) 


则 




Pa 口 —六 o c ’-”)) 

1 d 


WBSSBr 


A 




o da 


+ € t ( x ^ — x ( ^) d<x 


IT 


S (々 


,(v— 


1 


”)〜 \ hx k (^ 


+ a(x {p > — x^ 1 } ))da > 


因此 


n n n 


iu - acx^ac^^w <s S Si^n^i 


< S S S - j p 

J —— -1 —— -i r _^. -i b i J 0 




=^1 fe — 1 . 


+ Cf ( 尤 


Cp ) 


））am 


< m II4 U) - ^r 1> IIIU(^- 1) )~ 1 [l 


( 10 ) 


这里用了 


2 l A；t U ^ \ S M 


命 


* 144 » 


PDF 文件使用 11 pdf Fact or y 11 试用版 本创建 www. f i nepr i nt. com, cn 
















m 


1 — A (x^) A = B 〈 11 ) 

则 

acyo )， 1 = aCx^-^-ki ― b)-\ 

由定埋的假定及上章 S 1 引理 1 得到所要的氣论. " 

定理3 如果 ' 

U f|^ (u - x^\\\\A(x^r i \\ <«a< - 

2 

<a 

及 

\ x iy) ― |j > 7 ， 



证明方法和第五章的定理3完全一样. 

在计算机上使用时，不必验证（5)，（7)， （8) 是否适合， 
而直接应用疊代公式（3)，如果 ； 

jjW + i ) — < I 卜⑴一 太 ( 叫||， 

即可大胆地做下去，不然，停机换其他方法缩小范围. 

我们提请读者注意第 五章& 5的方法. 

■ ■ 


S 7* 求重心 

在第三章中我 n 曾提到求重心的问题，萵维的重心必领 
应用多重积分，而多重积分的计算在计算数学中又是一件麻 
烦的事，一般讲来我们需要计算出/ + i 个 < 重积分 

F = j - - * J - • - dx ” Mi $ ■ . . j -…， 

R R 

u s == j …! * 必 “ 
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谜里 R 是一个多面体，也就是造合 P 个不等式 





^ (j — 1 5 2, 



的点的集合，计箅多重积分时 5 我们用 Mont - Cdo 方法，对 
这一问0有稍为简化的建议，假定 i ? 在一个“方体” 

K x \ x t \<A 〔；二 1，2,… 〆 ） 

内，利用随机数表(:或随机数发生器)在获得 f 个数（匕，■〃， 

■ 

6) 后，看它是否在 a 中.如果不在中就作罢；如果在 ie 中 
則写下 

U ， .，.，§/; ⑴ 

再获得… s ❻后，如果不在 S 中则作罢，如果在 
只中 ，则在 （1) 上加上 

H 妒，…，❿； 


'得 


2 ,巧卜 6 

再做下去，得如果不在 B 中，则作罢，如果在及 
中，则加成 

3，# = # +妒， ... W =^ + 

一般是 


〜 m 


( TO -1】 


=3 


巧产 -2)+| 产 —1)， . • • ， ^C^.^1) = 


■ 


当 TO 充分大吋，我们拿 


■fra —1) 


m 


Vs 


CfflZ —1} 




m 


来作为重心的坐标.何谓充分大？当然可以有槪率建议，我 
们的简单建 议是： 如栗再做一次和上一次所得的数值差不多 

就停止（或只要求和差不多就停止也司以，或 

热 rw + i 

要求 
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"+ 亡 > 

n 

差不多也可以). 

这方法的概率误差是极好的，因此可以认为已经符合我 
们的要求.也许有 人说： 概率总是概率，其中包括了万一的 
失误，但是，这一点小险是要冒的. 


-h 


#■ » 


Cm—1 》 


和 S 


严+ 


m 


尤 

I 

0 ， 147 1 
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第三部分附 录 


附录 _ 度置问题 

令中00是一个不降的函数，而且 

屮 ( 0 ) = 0 ，屮⑴ = 1 , 

则当 X 过 [CK 1 ] 时 0 =屮 0 )也从 0 到 1 , 

在变数 r 上用黄金分割法逐步在 

〜〜 … 

做试验时，等价于 Y 在 

如 A, y 3 , … 

处做试验，而 

由中值公式可知 

y • — y«-i ^ <p f (SK^n - 丈^)， 

§在太 m 之间，如果 < pXi ) c 1 ,则用 y 做变数比用 x 
好，不然则不好， 由于 〆 00的平均值 

= 中（ 1 ) 一 qt)( 0 ) = 1 + 

及 tp ' O ) > 0 ,所以 rp ' OO 有时大于 1 ，也有时小于 1 * 

这说叨了一个问题，如果函数/ 00 的最大值在则取 
< P J (^ c ) < 1 的)/ ~ ( pOO 作为测度为好，但由于我们不知道 
最大值的所在 j 因而无法取合适的 vO ). 也就是任何一种测 
度用黄金分割法〖有些/ 0 )上算，柯些 fOO 不上算，也 
就是用来回调试法，总会碰得巧上算，碰得不好不上算，但上 
算的可能性极小,因此在一般情况下，我们就用原来普通的测 

* 14 S * 
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度就行了. 

多变数的情形测度的变化影响更大，而且有相互关系，但 
在实际问题中，一般地用常用的标准度量单位就行了，不要乱 
改度置单位，但不排序有理论根据，而选取合适测度，例如在 
知道每种变化与几何级数有关时不访用对数刻度. 
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附录二与后道工艺过程无关的优选法 


在实际中曾经遇到过以下的多因素优选法：要求在区间 
( 0 ， A ) 中选取 4 个点 

a b 

使目标函数 b ) 取最大值，但客观上这问题 
有以下的特 性：在 a < x t - 中 ，〜也是使某一函数 

/ fC ^3 J — 15 最人值 的点， /‘(’ I ， * * * ? ^ V - 1 ? ^ 5 

与 A 以后的点的选择无关+ 

因为 I 是使〃， r a ) 取最大值的点^也就是 a 适合 
于 

足 i ( A ， ^ 冬 AOi ， *?， ^)^0； 

dx 1 

而 JT t ， ^也是使 / sC^u X 2 \ x 3 ) 取最大值的点，也就是适合 

于 

心 ( Uj ， Xi ) = — AO. ， tf ， ^)^0； 


因此变为求如下形式的联立方程的解的 问题： 

x 3 ) =^ 0 , 

< . . ⑴ 

g 卜 i(l ， - * *， \) = o J 

L g^( x ^ * h) = 0* 

先取一个试点 h = 4”，由 （ I ) 的第一式解得 h 由第 

二式解得 h •… 3 由第^ — 1个式解得4以 
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W 1 ， 代人 o ) 的最卮一式得 

™ …， 

如杲 h = 0,则 ■ • •，就可能是解了3不然不妨假定 
心 > o s 再取4 =十 3 同法得 

h 2 = gk(^?\ …， 4 °) 

r a> _i_ y ay 

如果 h =0, 则得解答，如果心<0,则取4= 1 十文 1 进 

2 

行计算，这样用对分法很快就能得出 （1) 的解的近似值.如果 
心>0 5 则比较屺乂哪个大，例如则在…产 
的相反方向选4 5) 就可以找到问题的近似解了， 


• m ■ 
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附录三 一致 分布点寻优法 

从本书一开始就假定了所考虑的函数是单峰的.在一般 
工业生产过程中及在理论上较明确的科学实验中，已经在生 
产中的工艺中，这2假定是可以保证的，这是我们所介绍的优 
选法能很快地取得大量成果的道理.但不能排斥，在有些认识 
不足的探索性的科研问题上会出现非单峰的现象，我们就遇 
到过“液晶”的配方问题.如果照“穷举法”或“选优法”做，工 
作量太大，几辈子都做不完，而且问题不是单峰的，但与一些 
化学工作者在一起灵活地混合运用了本书所介绍的方法，不 
但很快地找到了最优点，而且找到了一条“山脊' 

在实践经验少 ，不会 抓主要矛盾的条件下，在理论认识不 
够，不曾获得精确范围的情况下，尝试新的科学研究课題，以 
下的方法，可起辅助作用. 

如果求々个变数 U ，• • ■ ， T 的函数 / O ) 的最大 
值，其区域是 

Si < < Vi (― 1 ， 2 ， … ，々 ）. 

在区域 D 中取一个一致分布的点列 

/’） （卜 1，2, 3，…乂 ⑴ 

在做了若干次试骑得出 

f (，) (卜 1 ， L …， iV ) 

后， 可得相当满意的函数值.再在其附近用本书所介绍的优 
选法. 

这一方法的优点在于容易发现复杂的"地形”（非单峰的 
情况），但缺点在于虽然比选优法好得多，但比起优选法来次 
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数还嫌太多， 

至于一致分布点列的构造，就不在此叙述了.请看华罗 
庚，王元《数论在近似分析中的应用 》( 科学出版社， 1978). ^ 
致分布点法，在求数值积分上已有好的成效，而求极值的问题 
可以看成为某一带参数的积分的极限，即 



PDF 文件使用 11 p df F ac 1 0 r y 11 试用版 本创建 www. f i nepr i nt . com, cn 








附录四求定正方阵的逆 


在第二部分第五章 H 中要用到求定正对称方阵之逆， 
同时定正方晖求逆法也是种类繁多的多因素优选法的根本依 
据，因此写这个附录，先总的指明一下，错综复杂的许多方法， 
实质上都是从解二次方程的“凑平方法"得来的. 

1- 求对称定正方阵的逆 



是一对称方阵 ， S = V ，其第 *1 个主对角线方阵是 


/ ^iL? * " ^ Hi 



又令 = •. ■，心，卜上乂的逆巧（-^广）可由以下的 

递归公式得之： 

5 f l = = 丄， 




一 1 
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( 1 ) 


(-P^ "㈣， d 

- 1 ** ■■ — ■- ^ ■■ — - - - 1 —— 

〜一 V ^ H^VU 

因此经々次迭代，得 



^1-1 


0 


0 S “ 一 1 MV ，: 

(即“凑平方法，，求逆得 


Sri = H ⑺ 


I — H ㈣ F ; 


0 




s “一 y^H^vU 
即得所证. 



-^ v^v 

- V^H U ^ i 



5 


附记求任意方阵 3 的逆也可用此 法：由 恒等式 

A~ l = AX^A T )- X ? 

化为对称方阵的求逆法，有些书对3 时, 

将 AXAA r y x 定义为 Z 的广义逆. 


2* 线性方程组求解 

在第二部分第五章经常遇到线性方程组 


M 55 _ 


PDF 文件使用 11 p df F ac1 0 r y 11 试用版 本创建 www. f i nepr i nt. com, cn 









xA =^b 

的求解问题.这里 j 是 〃阶对 称正定方阵，&是》维向量， * 
则是〃 维解向量. 

射于这种方程组的求解 ; 建议采取以下的方法. 

第一步：先将系数矩阵 j 表成 

A - LDL T 

的形伏这里 L 为一个对角 线元索 全等于1的下三角阵，1/ 
为其转置，/?为一个对角矩阵. 

用 ~(1 < ^ < ») 表示矩阵 d ， L 中位于第 〖行， 
第1 列交点上的元素，自然有 

〜 — 攻"， 

hi ^ 1^ hi ^ 0 (1 < / < 0， 

又用 < 表示矩阵 D 的笫〖个对角线元素，于是根据公式 

hi - - S o^o - hk )(1 

|E -1 ' 

an — 、 hk 在 k)lik O < ^ S 

k^\ 

依次确定出 L 与 D 的 元素： 

■ ■釋 > * * * 

^ Sriy 纟 “ 3 ， * * * ? - 1 ， ^ iy 

HP * V- «-*-*• *_■ * ^ % 

^*1 ? ZflJ ， * * * tt—i J ^nm 

在实际计算中 3 对于一个固定的；由于 

(1 1 ) 

的反复运用，可以在每舁出 一个心 后乘以 A 而贮存起来，以 
便后面取闱. 
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易见这一部分的计算 ft (乘除法)的数量级是 W / 12 . 
第二步，在确定了 i 与 D 后，通过解方程组 

zL f ^ b 
JtD = JB： 
xL «= y 

就得出原方程组的解. 因为 L ， D , X /都是三角形阵或对角 

线阵:所以这些方程组都很易解,其乘除法计算量的数量级都 
不超过 oQn 1 )• 


3 . 求抛物体的极值 


令 


Z 


a 


私十 


■ II I k - 


2 



⑴ 


此处 P 是4维向量， S 是 k 行列的对称定正方阵、“凑平方法 53 
就是以下的 方法： 


z a — ~ pS~ l p + 丄 O 屮 pS-^Six H- pS^y 

2 2 

^ a _ -r * 〔3) 

这就是当且仅当 

x — — pS - 1 (4) 

时 3 = 取这最小值 a — 丄 

2 

如用 L (1), 则得解方程 （4) 的方法如下， 

命 P = ( Pi ，■ . .， Pfdy P U：f ^ C / 7 ^ …， Pi ) 定义 

=■ p a> " ⑴， ..■ 


X 


C/) 
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p u K - 



u—' iy(—Pwff 0 - 1 )， o 

S “一 

H —— ~ C， ~ n ^-7 十 上 -(— F — ， H 。 -” ， 1 > 


而就是 （4) 式的解，这就是求 （2) 式极值的方法， 

但这样计算置嫌大,我们可 以用： S * 的办法去求方程 

xS = -p ⑴ 


的解. 


附记有些书上讨论矛盾方程组 


xA = p, 

龙是免维向量， P 是 Z 维向量，4 = ^ < /的求解问 

题.这神方程组本无解，但硬以 

作为其广义解 ，应 用时请注意，这仅对已知 x 与 p 确有线性 
笑系，并且有观测误差时，才可能有用. 

如果 x 与 P 并无线性关系，而滥用此法，不但不能正确 
地反映客观现象：反而导致思想混乱.例如原来客观问题是 
在某一范围的内点有极值，如果由若干点的试验数据，而用求 


广义逆的方法(或线性迥归)来处理，由于线性函数在任一范 
围的内点是无极 值的， 很明显这只会导致错误的结论. 




4 
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附录五离散与连续 



1- 在第二部分第五章中所介绍的抛物体法是离散型的 
方法，如果我们所考虑的函数是易于求出其一、二阶偏微商 
的，则有以下的连缤型的方法. 

以太=(4^ ^ - , ^0表示々个变数的向量，假定函数 
fOO 有一、二阶偏微商，以 


▽/GO 


( 1 ^ 



表示 （ 的梯度向量，以 

、 dXi - i<k 

表示 f 的二阶微商矩阵(又祢陈），则由在一点 
«处的幂级数展式的二次以下项 


/( tf ) + ▽ f ( a )( J ： — * d )’ + 丄（，;一 a )' (1) 

2 

可以得出迭代法. （ 1 ) 式的极值可由“凑平方法”得之，如在 
附录四中所做的那样，若^^为定负矩阵，则 （ 1 ) 式就是 

f ( a ) + 丄 O —这 + — a + VfH^y 

2 


—夺 v / Hr 1 ( vO , < -丢 v / w / Xv /)% ( 2 ) 

2 2 

即 （1) 式当 

x ™ a — VjHf 1 (3) 

时取极大值.令 r = 是一个初始值，则用递推公式 

4T U+1) =- — Vfix^HyXx^) (n = 2 3 -*> 
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比第五章更容易证明这方法也是 log log 型的方法. 

2, Davidon - Fletcher - Powell 引进了如下的半离散半连续的 

方法，从初始值及出发，用以下的迭代公式 M 

得々十”， h a+u 

有了 ^ rt > ? H n 后，求单变数 A 的函数 

印⑴司 — 入 ( 4 ) 


的极大值，得 2 

= K . 然后取 



并令 




及 



H fl+L = // n A 


(▽/( 

vKx^yn^y ： 


(yM,y(yM 

y rJ^nV rt 

关于这个较复杂的收敛性问题，我们只得到 


X 


<0) 


|卜善 ⑷一 D ， 


⑴ 


㈤ 


⑺ 


这里，是 K 幻取极值的点， （7) 式仅仅比 

II ，) 一 ，州 !I ◊ < 1 

II ， 一 ：<叫! 〜 q 
稍好一些，而后者是属于 bg 型的， 

0这方法中要求梯度，如果用差分来代替，则和本书第 
二部分第二章所介绍的最陡上升法有同样的缺点. 

2 ) 的极值怎 样求？ 如果用单因素优选法 ，一 次就 

要用 bgbg 次试验，仅就这一点求法就不可能在数量级上好 


过抛物体法. 

3) 求梯度是连续性的， 而求化 是离散性的，如果9>(义) 
求极值是用解析方法，则自然也就假定了 / O ) 可以求二阶偏 
微商，因此 h 可以用 O u 〜 一/， H 来代替(相差不多），这 
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里 H 是 K#) 的二阶微商矩阵 (Hessian matrix). 

4) 不要因为他们 只能估 计出收敛速度(7)，而以为 （7) 就 
是本方法的评价，实质上，用本书的方法是可以估计出比 （?) 
更精密的收敛条件，但不再把时间花在不比抛物体法好得多 
的问题上，因而作为读者练习（提 示： 请考虑 Ik ⑼一 > + 叫1 
与||^>__^||的关系乂 

3 - 在1•中介绍了一种连续型的方法，这类方法无疑可以 


改进.例如我们还可用以下公式 


, u+ij ^ 


； df ( x ia} ) 


dxi 


2 




df ( x u Q 


df ( x w ) 


O = 1，2，…） 


Qx ^ dxpdx ^ 

进行迭代，这里上下有相同指标表示从 i 到丨求和， 


O 严 - 4"))， S = . - //(/(>("〕））一 1 ■ 

W …衫/ 

这个方法比 1. 的方法收敛得 快些， 但还是 loffbg 型的 方法， 
且计算量较大，与这方法相应的离散型及半离散半连续型的 
方法，由于既没有理论上的困难，也没有实际上的显著优越 
性，这里就不谈了. 
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附录六几何优选法 


1.光以最短的时间从一点到达另一点.两点之间以直 
线为最短，因此光从一点^到另一点 S 的行程为连 B 的 
直线. 

光由一点^经镜面上0点反射到另一点图中角《称 
为射人角，角？称为射出角.如果 a 今？，作 B 对镜面的对称 
点 OA^OB > AB \ AB f 交镜面于 O 、 则显然有 CM + 

OB > 0 f A 4- 因此由光行时间最短推导出射人角《 — 

定等于射出角 /?. 




光的折射也是这个原理*例如光在空气中以速度^进 
行，在点 C 进人水中以速度“达到 B 点，求最速的途径.取水 
平面上一点 O 作为原点,./， C 的坐标各为 （〜， 〜乂（6, 


6)，0,0)，因此由点3到点5的时间等于 


/0O 

ii 
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求 / U ) 的微商使之为0, 


a 


x 


夕\/ ( 


-* A *) 2 H 1 - 


by — X _ 

u%] Ch l — x^ 2 i>2 



因此得 


sin & 


f/ 


Sill 




即得折射定律. 

这个折射定律实际上就是光学设计的主要物理根据 .一 
个光学镜头（如望远镜，照相机镜头等）通常是由数片透镜组 
成(图 64). 来自镜头前物面上的光线，按照折射定律依次通 
过各镜面，在镜头后成像.但这个像不一定是原物体的一个严 
格的几何相似形，而是有一定的偏差(称为像差).所谓光学 
设计，就是在某些约定条件下，通过选择这些透镜的镜面曲 
率、厚度、间隔、折射率等参数，使所产生的像差尽可能小一 

些暑 





m 幻 


rV-f r=T+^ J-lT. 甘 n_l，. txj. t 

. j* 
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i ~ I i =2 1 1 - 


而距离 f 等于 


n 


fj 


X 




tanoc 


^ iP^i 


y ^_ 

T^I '/• i 一 〆 


tr 


这里 p ■，内 为第 / 层的波的入射角， 

如果地层为连续介质，且速度为深度的函数，则时距之间 
有如下关系 

t ^ 2 [ 2 和― : ^ = 2 ( 2 

^ v{z 、 4' 一 pv(o’ J Vi — ? vu)’ 

这里 s 为反射层的深度. 

之集散点问题.在 （ t , y f ) 点有〜吨原料 0 -1, 
2, •…， /) 要运到处加工，加工后又要各运卜吨到 
A ) 处 （/ = 1，\…，乂问选择哪一点 （6 »/) 运输的 
吨公里最少？吨公里等于 


— !) 2 + (:^ — ”) 2 


m 


S — i ) 2 + (q — » i ) 2 . 


要解联立方程 


a id — ^) 


0/ _ 十 

i n ■■ ^ ■ ■ F h . ■ 

dg sjixi — + G 〜 


V) 


2 


m 


s 


^;(l — j ) 


=1 一 l) 2 4 - (t” 一 7j) 


2 




s 


“iiv — Vi ) 


v( A — 妾 ) 2 + (h — 巧 ) 1 


0 
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nt 








0*i — f) 2 + (乂 ， ? 一 ??) 


i 



中的5, >7 的数值并不容易.1960年万哲先同志提出了以下 
的模拟法. (注 意: 当/ _ m — 1时，这间题就已比上节光 
的折射问題难些了 .） 


先在一抉光滑木板上绘制 f 个发点和 m 个收点的分布 
图，然后在这些点处各凿一光滑小洞.在板面上置一光滑金 
属小圆环，小圆环上系/ 根绳子，每条绳子各穿过一个小 

洞，通过发点的绳子系上质量为〜的重物，通过收点 
(力，的绳子系上质量为&的重物.这样当小圆环达到平 
衡位置时，平衡&置就是所找的7?乂 
我们引人平面向量 


«x 




X 


V 一 ？) 2 + {yi — 7}} 1 


V 


V 




y -- - "»■ -- — - ■■ -- ■ - 

V (A — f) 2 4- ( y { — 7 i ) 


0 *i 一 f) 2 + (a — jj) 


2 




Vi 


_ f) 2 + ( 〜 一 TJ) 


则上面的联立方程就是 



占 A 



+ 2 = 

J - 1 


这正是/+%个力的平衡 方程， 也就是上述模拟法的理论根 

据. 


3. 邮路问题.一个邮递员每次上班，要走遍他负责送信 
的每条路，然后回到邮局.问应该怎样走才能使他所走的路 
程最短？这就是邮路问题. 
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首先考虑这样一个问题.如下面这样 一 个 道路图（每 
条线上的数字表示距离，小圆圈中的数字表示交叉路口的编 
号）.我们从一点出发，有没有可能沿每条路都走过一遍（不 
许重复也不许遗漏），然后回到出发点？ 



我们抝交叉路口分为两类，偁数条路的交叉点称为偶点， 
奇数条路的交叉点称为奇点.例如上图中，①，@，⑦是偶 
点，②，⑧， ®, ⑥是奇点.易见，当一个图中只要有一个 
奇点时，就不可能那样走一遍了，因为对一个路口来说，有走 
进去的一条路，必有走出来的一条路.反之也不难证明，当道 
路图上没有奇点时便可以那样走一遍了.如果起点和终点可 
以不同，则可以允许有二个奇点+这就是所谓一笔画问题. 

回到邮路问题，如果道路图上没有奇点，问题就已经解决 
了.如果有奇点，我们总可以通过在某一些路上再重复走一 
遍的办法来消灭奇点，这相当于在图上再添上一些路.上图 
中， 如果在@与③之间， ® 与®之间各添上一条路（即在 
②一③，⑤一⑥上走两次），奇点就没有了，于是就可以一笔 
画了，但我们还要求所走的路程最短，因此问题归结 为：如 
何选择要重复的路程，既能消灭奇点，又能使重复的路程最 
短？ 

这个问题已被管梅谷同志解决.他证 明了： 只要在每个 
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圈上董复的路程不超过整个圈长的一半，这时就可以得到路 
程最短的走法了， 

例如在上页图中，如果重复@一@，⑤ 一⑧ 两条路，可以 
消灭奇点，但在圈③一⑧一④一⑤一 ⑥一③ 上重复路程2 + 
3 =5 已超过整个圈长之半 ( 4 . 5 ), 所以它不是路程最短的走 
法，但我们可以改为重复② 一®， @一④，④一⑤，这样既 
消灭了奇点，又使得每个圈上重复路程都不超过整个圈长之 
半，于是得到这样一个路程最短的 走法： 

綜合上述，得到谊样一个解决邮路问题的歩骤：首先把 
道路图上所有的奇点找出来，然后选择一些需要重复的路，使 
奇点消灭.再在每个撣上检査重复路程是否超过整个圏长之 
半，如有超过的情况，则在这个圈上，把原来不重复的路定为 
重复的路，而原来重复的路定为不重复的路.这样逐步调整， 
当所有圈上重复路程都不超过整个圈长之半时，就达到路程 
最短了.最后将整个图形一笔画出. 

一般来说一个图上所具有的圏数是很多的：是否必要每 
个圈都险査？如上图共有六个圏，实际上只要捡查三个圈就 

够了. 一般，如一个图有 f 条路，》个点，则需要检查的圈数 
为 j + 1. 

-在画交通图时还要注意一个问题，一条很宽的马路，实际 
上是需要在马路的两侧各走一次，这时在交通图上要画两条 
路，而不 是画一 条路. 

这里提出的问题，不只是邮連员会碰到，在其他一些场合 
也会碰到，如马路上扫地、喷水的车子所走的路线， 

4. 蜂房问题.我们研究与蜂房结构有关的一个数学问 
题，从正面看，蜂房是一个正六边形，但它并不是一个六棱 
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柱，它的底部是由三个相同的菱形拼成的（图 67). 菱肜的钝 
角为109°28%锐角力70^2^说得更具体些，拿一支六稜柱 
的铅笔，未削之前，铅笔一端的形状是 / fBCDEF 正六角形 
(图 6 S ), 通过—刀切下一角，把三角形 ABC 搬置 AOC 
处， CE 切如此同样两刀，所堆成的形伏就是一个蜂 
房. 一 个蜂巢就是两排这样的蜂房，底部和底部相接 而成， 




建造一个蜂房所消耗的材料是与表面积成正比的 * 我们 
问： 怎样切出来使所拼成三个菱形做底的六棱柱的表面积最 



小. 

假定六棱柱边长是1，这吋 

AC - VT , 把图67的尖顶六棱 
柱表面分成六份，把其中之一摊平 
下来，得出图69的形状，从宽为 
1的长方形切去一角，切割处成边 

以处为腰 ， ^为高作等腰 

2 


三角形，假定被切去的三角形高为 h 则 

ap^ 
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pp 


,=2 V 


x 


2 



1 


V 1 4 - 4 欠 2 , 


^ APP f 的面积为 


VIVT 


4 


4 - < 


所以问題变为求; r ， 使 


fOO 




2 


x 



3 


4 



-h 


2 


最小 • 

我们用一个初等方法来求 fOO 的极小值.令 2 jt 
^0 >0), 则因算术平均 > 几何平均，可见 





3 




t=i 


V 3 — 1 



3 


w- 


4 


4l 


>2 


/(V 3 — 1)(V’ 3 + 1) _ _ 1 

V 4 3—7 


而旦仅当 



3 一 1 



3+1 


4 


斗 


4/ 


时取等号，这时 


v ； 


3 


x 


V 

飞 


，，而 


2 


2 VT 


r ■ ■ ■ 

2 V 2 


4(1 + \/ 3 ) f V 8 


所以 x 



时 wo ) 取极小值 j 

V 8 
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令 = r ， 由余弦公式得到 
2(1 4- jc 2 )cosr = 2(1 + x 1 ') — (1 4- Ax 1 ) = 1 — 2x 2 ^ 

UU = | 代 A 得到 

V 8 

1 

CDS Y 5=31 - - 

2(1 

因此得到 r = 70°32'，与蜂房的实际角度相一致.我们进一 
步问： 蜂房是否也具有“在体积一定时表面积最小”的特性? 
下面将 证明： 一个形如蜂房的尖顶六棱柱，体积 P 为一 



定时，表面积(不算底面）最小值为3^7 F 2 、 而且当六角形 
边诠是 為 高度是 ^ P 时取最小值. 

以边长“的正六角形为底，以 * 为高的六棱柱，其六个顶 

■ ■ - 

点顺次沟如上述，过谷（或 D 或 F ) 棱处(假 

*v 8 


定$ > 



心及〆， c (或 C ， E 歲 E ， A) 作一平面 5 切下 


三个四面体，反过来堆在顶上，就得到一个以三个菱形做底的 
尖顶六棱柱， 易觅它 的体积 P 和表面积^为： 

V = hlj^a 2 b, 

2 

S 6 ci (jf -h 


我们有 







= 3VT V m , 
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而且仅当 


2V 


V 3 ^ 

时取最小值，这时 

b s = : — 



3 





2V 


3 V 3 


V 


j /切 



■ 


尖顶六棱柱高度为 








2^ 


171/3 


*■ 


即上述结论得证. 


通过实测，蜂房的大小与上例并不一致经过实测， 


0. 3 5 cm ， 高为 0.70 cm , 而按上面的结论，高应是占 


4 






^0.38 cm ), 所以蜂房并不具有“在 


体积一定时表面积最小”的特性.而蜜蜂是根据自己身村设 
计，用材最少的窠. 




_ 
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附录七多目标问题 


在实际中经常会遇到多目标的问题，也就是我们要优选 
一 些参数5使几个目标函数同时都达到较好的水平.这里我 
们建议这样一个方法，在这些指标中，根据当前的情况，选择 
一个优先要考虑的指标作为我们的目标，而其他的指标希望 
它不要趄出某个界隈，也就是引人一些不等式的约束条件，这 
样，问题就化成单目标的问题了. 

有些实际问题可以把多目标问题化为单目标问题.例如, 
用一个关于经济价值的指标把原來关于劳力、各种衬料消耗 
等的指标合为 ■一 个. 
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附录八重复试验 




在得出较好的方案后，按照这个方案固定所有因素反复 
做试验有时是必要的，但可能会得出不完全相同的数据，对这 
些数据怎样处理？如果做的试验次数不太多，例如只做几个， 
十几个或几十个，可用以下的 方法： 

把试验数据按次序排列起来 

得 结论： 就做几次试验的结果来说，100%落在范围 
之间，落在间的可能性是等 
等， 

如果做的试验次数实在多，可用以下的方差分析 法：先 
求平均数 

I = ™ C^i + * * * + 

n 

及方差 

■ s 2 =- 丄 7 E (元 一 a ) 2 ， 


而结论是：在 b 


s 


邱讳；在 U 





% 4 - 


% 


S 



之间的可能性约为 


n 


2S 


n 



之间的可能性约为等 

n 


等.（附记为什么沪的分母是《 — 1，因其分子为 


■ra 




•■J --7 _ 

? 一 nx x = 




(乂 + 


f ip _ 


n 


+ O 2 , 
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这是 A ， 的二次型，它的方阵 

/ — 1(1，…，1)'(1，…， 1) 
n 

的秩显然为《 — 1，也就是实际上 f 的分子不是《个数的平 
方和而是《 — 1个数的平方和). 

如果得出两种较好的方案，各自固定所有因素后反复做 
试验，可能会各自得出不完全相同的试验数据，怎样比较这两 
种方案的优劣？如果做的试验次数不太多，则可将所有数据 
按大小次序排列起来，按本书第一部分第三章 S 9所介绍的 
办法来判断.如果做的试验次数实在多，可用以下的方差分 
析法： 若第一种方案做了 A 次试验，得数椐为〜， 

第二种方案做了 h 次试验，得数据为 • • ■ , 分别求 
平均数 

+ - 1 

5=1 — Xj ^ ~ x i 

n i t = \ / i 

及方差 

^ —- 2 d — 文 ") 2 ， 幻在 —— 2 d — %) 2 . 
—— i 薄 3 一 1 i==1 

若 \ x 3 - x { \ > 1,645 /^1 + ^, 则约有90%的可能性，这 

V 

两种方案所得的结果是有差 异的; 若 

\xi — ij > 1.96 ^ J ■ ， 

则约有95%的可能性，这两神方案所得的结果是有差异的； 
等等. 
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附录九 0 - 1 变元法 




1. 在优选法形成一数学分支之前，曾有一些 方法， 科学 
根据不太足，先验地假定某一参数仅在若干个水平中选取，最 
简单的是假定在两个水平（或三个水平）中选取，非此即彼.他 
们所甩的数学工具，与其说用了连续变元，不如说是周了 0-1 
变元，即取不取变元，以两个水平为例，温度或取300^或取 

我们不妨说取300^时这变元的值是1，不取 300 ^C 

时，这变元的值是0,换变数就变为 tM 奎元 

600 300 

q 了.下面 3 我们简略介绍一下5并且从原则上更提高一些. 
一 个参数$限定在”个水平 

§1， li ，• ■ ■ 5 

中 选取其中之一 ，可以 用《个 bl 变元 h ， 来代表，送 
捏;一 •，:只取值 Q 或1 3 而且 


即《个水平中每次只能取一个. 

如果某一目标函数 T , 依赖于4个参数 
而泸> 有〜个水平 


泞)， •. 




其对应的 0-1 变元 …， 假定其自标函数: T 是最简 
单的一次模型(这种假定的可靠性极 小）， 

^ = (o 

i = \ j -1 
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2>5 0 =1， # = o 或 i. 
/ "1 


用 








代人 （ 1 ) 则得 




⑴ 


+± sc ^ - o 尸. ⑶ 

i == 1 i = I i = i 

这式子共有〜— <个独立的 o - i 变元 j ^， 但必 
须适合 

， ( 4 ) 
/ = 1 

在 k 个参数中，每一参数取一个水平的情况做一次试验， 
也就是取一个 g -1 变元组得出一个 r 来，做了 n 次试验 5 
得出况个 K 

S 岭)一 o4? (1 <i 1 <; — 1) (5) 

i - 1 

为待定系数的一次联立方程组，待定系数总共有〜+ ■•• + 
- 々+ 1个，因而做 

N =、冬 4^«是——々4^1 (6) 

次试验，就可以从解一次联立方程中得出 （ 5) 的数值.但请 
注意，我们所做的试验，要使所 ffi 到的行列式不等于 0 ,定出 
(5) 后，就用以下的方法判断参数 £0 该取哪一个水平.这方 
法是，如果对一个参数 

0^—4? (; — 1 , …，〜 一 1) (7) 

都是负的，则取1,即0的第〜个水平如果不 
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然，而 （7) 中以/ = ； o (0 为最大，则取 =】，即取水平 
^ ^ 因此按模型，^在 

f C 1 )- 找;0,…，# = §軋 

处取最大值.因而他们就得出结论，这一点是目标函数: r 的 
最优点了.当然模型是很不可靠的，因为线性模型是没有最 

大值的，有些同志认识到这方法的不可靠性，祀这一点称为有 
参考性的最优点，也就是承认了这一方法不属于有收敛性保 

证的优选法的范围. 

2 . 为了更明确地说明问题.，我们以每个参数仅取两个水 

平为例.例如参数0只取两个水平在 P 或纟 P 之一.换变数 

一妒 一 W 
，杜) 一沪， 

则当取第一水平时，^ =1，当取第二水平时（即当不取 
第一水平时 h ^ = 0 . 如果假定目标函数原来是 

7" =砷 ~h 芝] ct ^ a \ 

i = l 

也就完全等同于目标函数为 

k 

丁 = fio + Z 

f ■= ! 

的 o - i 变元^的问题了. 

由于基本假定易于失真，我们就不再烦琐地叙述多种多 
样选取 … ^彳）做试验的方法.更不谈如何用统计 
方法来分析其试验结果了，要之，本方法不但不能保证得到 
最优点，甚至于不能 保证自 己先验预定的 V 、％个点中 
的最优点. 

至子多做几次试验，即上节中取 V > + ■ ■ ■ - t - Wk 一 

K - hU 用附录四中所提到的求矩阵的广义逆法来解矛盾方 
程，那就更不对头了 I 
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例 1 取（: r ” * 4 S ^) = #5 ** 这里 

A 是一个4维向量，除第；分量是1外，其它诸支量都等于 
0. 所对应的4+ 1方阵是 


1 

0 


1 

1 

.p 1 

1 

* 

0 

1 … 0 

h 

4 

h 

0 

■i 

■ 

_ 

W 

* 

譬 

警 


1 

0 



这一方法对应于从一点出发，向每一方向都试探一步而后定 
上升的方向 ，很 明显这并不比一步一升高的“瞎子爬山法 〃更 
好些. 

例2考虑某一紀方由三种化学药品所配成，每种药品 
昕取的用量有三种水平，混合后加温，温度也是有三个水平， 
这是一个四个参数的问题，三个参数表示用量 ，一 个参数表示 
温度3因此要求做3 X 4-4+ 1 =9次试验，就能定出其线 
性目标函数，但这一数学结构易于失真，祈以也就不再讨论这 
批试验点的安排问题，更不考虑不分主次，把所有的因素善同 
地安排实验的方法了. 
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